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Abstract: We develop a general approach to the calculation of target mass and finite 
t = {p 1 — p) 2 corrections in hard processes which can be studied in the framework of the 
operator product expansion and involve momentum transfer from the initial to the final 
hadron state. Such corrections, which are usually referred to as kinematic, can be defined 
as contributions of operators of all twists that can be reduced to total derivatives of the 
leading twist operators. As the principal result, we provide a set of projection operators 
that pick up the "kinematic" part of an arbitrary flavor-nonsinglet twist-four operator in 
QCD. A complete expression is derived for the time-ordered product of two electromagnetic 
currents that includes all kinematic corrections to twist-four accuracy. The results are 
immediately applicable to the studies of deeply-virtual Compton scattering, transition 
7* — > Mj form factors and related processes. As a byproduct of this study, we find a series 
of "genuine" twist-four flavor-nonsinglet quark-antiquark-gluon operators which have the 
same anomalous dimensions as the leading twist quark-antiquark operators. 
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1 Introduction 



As well known, target mass corrections ~ (m 2 /Q 2 ) k to the structure functions in deep- 
inelastic lepton-nucleon scattering (DIS) can be calculated exactly in terms of the leading- 
twist parton densities. On a technical level, these corrections arise because of the sub- 
tractions that are needed to form traceless operators, as demonstrated in the pioneering 
paper by Nachtmann [1]. The Nachtmann power corrections have been studied in much 
detail (see e.g. [2]) and are routinely taken into account in the analysis of experimental 
data. Another classical result is due to Wandzura and Wilczek [3] who have shown that 
the twist-three structure function g2{x,Q 2 ) for massive targets with spin 1/2 receives a 
contribution related to the leading-twist structure function gi(x,Q 2 ). This contrbution 
appears to be numerically dominant for the nucleon, see e.g. [4] for a recent analysis. The 
Wandzura- Wilczek relation follows from Lorentz invariance and can be understood as spin 
rotation in the target rest frame [5, 6] . 

The target mass effects are usually referred to as kinematic corrections since they can 
be expressed in terms of leading-twist parton distributions and do not involve "genuine" 
nonperturbative effects due to quark-gluon correlations. In this paper we develop an ap- 
proach to the calculation of kinematic corrections to a larger class of hard processes that 
can be studied in the framework of the operator product expansion (OPE) and at the same 
time involve momentum transfer from the initial to the final hadron state. This, more 
general, kinematics is relevant for exclusive scattering processes which involve off-forward 
matrix elements of the type 

tf\<D\p) t p 2 =p> 2 = m 2 , t = (p'-p) 2 

and form- factors at large momentum transfers (and weak decays of heavy mesons/baryons) 
that involve hadron-to-vacuum matrix elements 

(0\O\p), p 2 = m 2 , 

e.g. 7* — > rpf, B ->■ K*lv t etc. 

The forthcoming studies of hard exclusive scattering processes at the future 12 GeV 
facility at Jefferson Lab provide the main motivation for our study. It is generally accepted 
that such reactions allow one to access a three-dimensional picture of the proton in longi- 
tudinal and transverse plane [7], encoded in generalized parton distributions (GPDs) [8, 9]. 
One of the principal reactions in this context is Compton scattering with one real and one 
highly-virtual photon (DVCS) which has received a lot of attention. The QCD descrip- 
tion of DVCS is based on the OPE of the time-ordered product of two electromagnetic 
currents where the GPDs appear as operator matrix elements and the coefficient functions 
can be calculated perturbatively. In order to access the transverse proton structure one 
is interested in the dependence of the amplitude on the momentum transfer to the target 
t = {p 1 — p) 2 in a sufficiently broad range. Since the available photon virtualities Q 2 are 
not very large, corrections of the type ~ t/Q 2 which are formally twist-four effects, can 
have significant impact and should be taken into account. 
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The problem is well known and its importance for phenomenology has been acknowl- 
edged by many authors [9-19]. At first sight, such corrections are seemingly disconnected 
from nonperturbative effects (e.g. one may consider a theoretical limit Aq CD <C t <C Q 2 )- 
Yet their separation from generic twist-four corrections ~ Aq CD /Q 2 proves to be surpris- 
ingly difficult. 

The kinematic corrections ~ t/Q 2 and ~ m 2 /Q 2 induced by the subtraction of traces 
in the leading-twist operators, which are a direct generalization of the Nachtmann's cor- 
rections in DIS [1] were considered already in Refs. [10, 15-19]. These results are, however, 
incomplete, because ~ t/Q 2 corrections in hard exclusive processes (and for spin-1/2 tar- 
gets also ~ m?/Q 2 corrections) also receive contributions from higher-twist operators that 
can be reduced to total derivatives of the twist-two ones. Indeed, let Pl ,,, Pn be a multi- 
plicatively renormalizable (conformal) local twist-two operator, symmetrized and traceless 
over all indices. Then the operators 

are, on the one hand, twist-four, and on the other hand their matrix elements are obvi- 
ously given by the reduced matrix elements of the twist-two operators, times the momen- 
tum transfer squared (and for spin-1/2 hadrons also target mass corrections). Thus, the 
contributions of both operators in Eq. (1.1) must be taken into account. Moreover, as 
we will see below, the distinction between the kinematic corrections due to contributions 
of leading-twist and higher-twist operators is not Lorentz invariant and has no physical 
meaning; they must always be summed up. 

Again, at first sight, taking into account the operators that are shown schematically 
in Eq. (1.1) should not be difficult, at least to the leading order in strong coupling. The 
problem arises because 02 has very peculiar properties that can be traced to the Ferrara- 
Grillo-Parisi-Gatto (FGPG) theorem [20]): divergence of a conformal operator vanishes 
in a free theory. As a consequence, using QCD equations of motion (EOM) the operator 
02 can be expressed as a sum of contributions of quark-antiquark-gluon operators. The 
simplest example of such relation is [21-24] 

WO^, = 2iqgF v ^q , (1.2) 

-H- 

where O pv = (1/2) [qj^D^q + (/i -H- v)] is the quark part of the energy-momentum tensor. 
For simplicity we consider massless quarks. The operator on the r.h.s. of Eq. (1.2) involves 
the gluon field strength and, naively, one would expect hadronic matrix elements of this 
operator to be of the order of Aq CD , which is wrong. More complicated examples involving 
leading-twist operators with two derivatives can be found in [25, 26], e.g. 

\^O mfi = -12iqY {gF p pD a - D a gF p p + (a <-> /?)} q - 4d p q(^gF ap + j a gFp p )j 5 q 

8 ~ 8 ~ 28 ~ 

- 3 di3qj a gF aa ^q - - d a qj a gF a/3 j 5 q + y g a pd p qj a gF ap q , (1.3) 



where 

"15 3 

O pa (3 = Sym^ —q-i p D a Df3 q - -d a d p q-i p q 



traces 
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and Sym^ a p stands for symmetrization in the Lorentz indices. The general structure of 
such relations is, schematically 

(dO) N = Y,ai N) G N ,k, (1.4) 
k 

where G^k are twist-four quark-antiquark-gluon (and more complicated) operators and 
the numerical coefficients. The subscript N stands for the number of derivatives 
in On and the summation goes over all contributing operators, with and without total 
derivatives (so that k is a certain multi-index). The same operators, G/vfc, also appear in 
the OPE for the current product of interest at the twist-four level: 

T{j(x)j(0)} t=A = CN,k(*) G N ,k ■ (1-5) 

N,k 

A separation of "kinematic" and "dynamical" contributions implies rewriting this expan- 
sion in such a way that the contribution of a particular combination appearing in (1.4) 
is separated from the remaining twist-four contributions. The "kinematic" approximation 
would correspond to taking into account this term only, and neglecting contributions of 
"genuine" quark-gluon operators. The main problem is to make this separation consistently, 
such that "genuine" quark-gluon operators are in some well-defined sense "orthogonal" to 
the kinematic contribution. 

The guiding principle is that the separation of kinematic effects is only meaningful if 
they have autonomous scale dependence. The twist-four operators in Eq. (1.1) are special 
in that they obviously have the same anomalous dimensions as their "parent" twist-two 
operators. Hence we can reformulate the problem in the following way: 

All existing twist-four operators with the same quantum numbers and of the same 
dimension mix with each other and satisfy a renormalization group (RG) equation which 
can be solved, at least in principle. Note that we have to include operators with total 
derivatives, so that this is a large matrix equation, in general. Let Qn,u be the set of 
multiplicatively renormalizable twist-four operators which we can express in terms of linear 
combinations of the original operators, i.e. 

GN,k = ^2 ^k!k' G N,V • (1-6) 
k' 

Relation (1.4) tells us that one of the solutions of the RG equation is known without 
calculation. Namely, there exists a twist-four operator with the anomalous dimension 
equal to the anomalous dimension of the leading twist operator, and Eq. (1.4) gives the 
corresponding eigenvector. (For simplicity we ignore the contributions of 8 2 On operators 
in this discussion; they do not pose a "problem" and can be taken into account relatively 
simply.) 

Let us assume that this special solution corresponds to k = 0, so that Q^,k=o = {90)n 
and ^1=0 k' = flfc ' • I nver ti n g the matrix of coefficients V'i^I' we can expand an arbitrary 
twist-four operator in terms of the multiplicatively renormalizable ones 

°N,k = ( d °)N + E A*' GN,k> ■ (1-7) 
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Inserting this expansion in Eq. (1.5) one obtains 

T{j(x)j(0)Y^ = £ c N , k (x)^ {dO) N + . . . , (1.8) 

N,k 

where the term in (80)n defines the kinematic correction which we are looking for, and the 
ellipses stand for the contributions of "genuine" twist-four operators. The problem with 
this (formal) solution is that finding the coefficients tp^J in general requires knowledge of 

the full matrix ip^^t, in other words the explicit solution of the twist-four RG equations, 
which is not available. 

In this work we show how to overcome this difficulty. Our main result is that the 
explicit solution is not needed; it can be avoided by making use of the symmetries of the 
RG equation that are rooted in conformal symmetry of the QCD Lagrangian [27] . 

Twist-four operators in QCD can be divided in two classes: quasipartonic [28], that 
only involve "plus" components of the fields, and non-quasipartonic which also include 
"minus" light-cone projections. Our starting observation is that quasipartonic operators 
are irrelevant for the present discussion since they have autonomous evolution (to the one- 
loop accuracy). Hence terms in {80)n do not appear in the re-expansion of quasipartonic 
operators in multiplicatively renormalizable operators, Eq. (1.7): the corresponding coeffi- 
cients vanish. As the result, the kinematic power correction ~ (0O)n is entirely due 
to contributions of non-quasipartonic operators. 

Renormalization of twist-four non-quasipartonic operators was studied systematically 
in [29, 30]. The main result is that with the appropriate choice of the operator basis, the 
corresponding RG equations can be written in terms of several SX(2)-invariant kernels. 
Using this technique, we are able to prove that the anomalous dimension matrix for non- 
quasipartonic operators is hermitian with respect to a certain scalar product, which implies 
that different eigenvectors are mutually orthogonal, i.e. 



E^«'«~4„ (1.9) 
k 

where /j,^ is the corresponding (nontrivial) measure. 

The knowledge of /ul is sufficient: using this orthogonality relation and the expression 
(1.4) for the relevant eigenvector, one obtains, for the non-quasipartonic operators 

where ||a^|| 2 = /4f^( a i ) 2 - Inserting this expression in (1.8) we end up with the 
desired separation of kinematic effects. 

The actual derivation proves to be rather involved. It is done using the two-component 
spinor formalism in intermediate steps and requires some specific techniques of the SL(2) 
representation theory. Some of our results were reported earlier in [31]. 

The presentation is organized as follows. Section 2 mainly serves to introduce the 
necessary formalism and notation. Following Refs. [29, 30] we define a complete basis of 
non-quasipartonic twist-four operators which correspond to Grjvfc i n Eqs. (1.4), (1.5) and 
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discuss symmetries of the corresponding RG equations. A new result in this section is the 
hermiticity of the anomalous dimension matrix with respect to a certain scalar product, 



Section 3 is the central one. Here we derive the expansion for divergence of the leading 
twist operators, (0O)n, in terms of non-quasipartonic quark-gluon operators in the con- 
formal basis (i.e. the precise version of Eq. (1.4)), and use hermiticity of the anomalous 
dimension matrix to reconstruct the coefficient in front of {80)n in the inverse expan- 
sion of an arbitrary quark-gluon operator in terms of the multiplicatively renormalizable 
twist-four operators, the analogue of Eq. (1.10). As a byproduct of this study, we find a 
series of "genuine" twist-four flavor-nonsinglet quark-antiquark-gluon operators which do 
not reduce to total derivatives and have the same anomalous dimensions as the leading 
twist quark-antiquark operators. For reader's convenience, we present a summary of the 
"kinematic projections" for arbitrary twist-four flavor-nonsinglet quark-gluon operators in 
Section 4. 

The next three Sections are devoted to the application of this formalism to the OPE 
for the time-ordered product of two electromagnetic currents, which is the most interesting 
case phenomenologically. In Section 5 we demonstrate how the twist expansion can be 
done conveniently in the spinor formalism, Section 6 contains the explicit construction of 
the kinematic contributions and in Section 7 we give the final result for the T-product in 
two different representations for the crucial twist-four part. 

Section 8 contains a short discussion of typical matrix elements that one encounters 
in physics applications. One issue that we want to emphasize here is that the distinc- 
tion between the kinematic corrections of Nachtmann's type, i.e. due to contributions of 
leading-twist [13-19], and of higher-twist operators calculated in this work is not invariant 
under translations along the line connecting the electromagnetic currents in the T-product. 
As we show on a simple example, translation invariance is only restored in the sum of all 
contributions. We also use this example to illustrate application of the leading-twist pro- 
jection operators and elucidate on the general structure of higher- twist matrix elements 
that contribute to hard exclusive scattering from a proton target, e.g. DVCS. 

The final Section 9 is reserved for a summary and conclusions. The paper also contains 
several Appendices where we explain some technical issues and details of the derivation. 

2 Conformal operator basis 
2.1 Spinor formalism 

The conformal symmetry which plays a crucial role in the further discussion becomes more 
transparent in the spinor formalism. The latter also presents a convenient framework for 
twist separation in comparison with the conventional vector formalism. Unfortunately, 
there are no standard notations (conventions) which vary from paper to paper. In this 
work we follow the notations adopted in Refs. [29, 30] which are summarized below. 

In the spinor formalism each covariant four-vector is mapped to a hermitian ma- 
trix x: 



cf. Eq. (1.9). 
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where = (1, a), = (1, — a) and a are the usual Pauli matrices. The Dirac matrices 
in the spinor representation read 

Here a^ v = | [7^,7"], 75 = Z7°7 1 7 2 7 3 and 

(O/ = \\^ V ~ , (^T/j = " frVfy • (2-2) 

The Dirac (quark) spinor q and the conjugated spinor q are written as 



X? 



q = bJ > ,$&), (2-3) 



where ip a , x 13 are two-component Weyl spinors, ipa = (ipa) , X a = (x a V ■ We accept the 
following rule for raising and lowering of spinor indices (cf. [32] ) 

u a = e a % , u a = u% Q , u a = , = e^vP , (2.4) 

where the antisymmetric Levi-Civita tensor e is defined as follows 

ei2 = e 12 = -e i2 = -e i2 = 1. (2.5) 

Note that for this definition ej = -e^ a = 5% and e"^ = -e^ 6 = 8$ and (e a/3 )* = e^ & . For 
arbitrary Weyl spinors we define the invariant products as 

(uv) = U a V a = —U a V a , (uv) = UaV a = —U a Va ■ (2-6) 

The gluon strength tensor F^ v and its dual F pv = ^e^ wp<T F pa can be decomposed as 

Here / a/ g and /.g are chiral and anti-chiral symmetric tensors, /* = /, which belong to 
(1,0) and (0, 1) representations of the Lorenz group, respectively. 
The explicit expression for the matrix x a a reads 



x + x 3 x± — ix 2 
X\ + 1x2 x - x 3 



x, L a» (2.1 



so that one immediately recognizes the two diagonal entries as the light-cone "plus" and 
"minus" components of a four- vector x^ [33], whereas the two off-diagonal entries corre- 
spond to the (complex) coordinates in the transverse plane. The spinor representation for 
four-vectors is in this sense similar to introduction of the familiar light-cone coordinates. 
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In order to retain the Lorentz covariance in the light-cone formalism one usually introduces 
two auxiliary light-like vectors 

n 2 = 0, ra 2 = 0, (nra)^O, 

so that arbitrary four-vector can be expanded as 

(nn) x,j, = x + h fl + x- n M + (nn) x± 

where x + = x^n^ and x_ = x^h^ . The same decomposition can be made quite elegantly 
in the spinor formalism by observing that a light-like vector can always be written as a 
product of spinors: 

"aa — Aq,Aq. j f^aa — Hal^a j (2-9) 

where X = X\ p, = fi'. The basis vectors in the plane transverse to n,h can be chosen as 
H a \a and X a fia- Thus 

(//A) (A/2) X a& = X ++ Hafia + X X a X~a + X_ + X a fia + X+- /Z a Ad , (2.10) 

where 

(2.11) 

so that x ++ and x correspond to the "plus" and "minus" coordinates, respectively, 

whereas x^ and x |_ are the two (holomorphic and anti-holomorphic) coordinates in the 

transverse plane. Note that in difference to Ref. [29] we keep the normalization of the 
auxiliary spinors A a , [i a and hence the product n^n M = (1/2) (//A) (A/2) arbitrary. This 
freedom will offer some advantages at the cost of slightly more cumbersome notation. 
Similarly, for quark and gluon field operators we define 1 

V>+ = A^a , X+ = A Q X« , /++ = \ a \ P faP , 

^+ = ^cX, X+ = Xa~X & , u + = f.fx«xP, 

■>!>-= H a il>«, ^_ = ^/2°, /+_ = AV/^ 5 (2-12) 

etc., so that, for example 

(A/z) Tp a = fj, a tp + - X a ip- (2.13) 

The "plus" field components correspond to "good" components that present independent 
degrees of freedom in the light-cone quantization [33]. For example, / ++ and /++ gluon 
operators can be expanded in terms of creation and annihilation operators for gluons with 
definite helicity, cf. [34]. Further details and useful identities involving the algebra of 
matrices can be found in [29, 30]. 

Finally we note that in a vector theory like QCD contributions of left-handed and 
right-handed quarks, ip a and x a , always enter in parallel. In the subsequent sections we 
mostly consider operators built of left-handed quarks; the remaining ones can be added in 
the final answers using general symmetry considerations. 



1 This "plus-minus" notation is, unfortunately, somewhat ambiguous because e.g. x++ = \ a x a a\ a but 
/++ = \ a \ f a p- A more consistent but less intuitive "lambda-mu" notation [29, 30] is to write ie++ = x^\, 
x+- = x X p,, /++ = fxx, /+- = hn etc. 
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2.2 Conformal symmetry and the SU (1,1) scalar product 

Despite being broken by quantum corrections, conformal symmetry of the QCD Lagrangian 
has important consequences and in particular appears to be very useful in the analysis of 
scale dependence of physical observables, see Ref. [27] for a review. In the applications to 
hard processes the so-called collinear subgroup SL(2, R) of the conformal group plays a 
special role. It corresponds to fractional-linear transformations of the coordinate along a 
certain light-like direction which is fixed by the kinematics of the process: 

i az + b / n iA\ 

x ll = zn fl , z -> z = — , (2.14) 

cz + a 

where a, b, c, d are real numbers, ad — be = 1. In particular we will be dealing with quark 
and gluon fields "living" on the light cone: 

q(z) = q(zn) , F^(z) = F^(zn) . (2.15) 

Representations of the SL(2,M) group T 3 are labeled by conformal spin j which can take 
integer and/or half-integer values. For a generic field (or function of the coordinate) the 
corresponding transformation law is 

^^ Ti ^ = jdw v &ii)- (216) 

The generators of the infinitesimal transformations corresponding to (2.16) take the form 

S+ = z 2 d z + 2jz, S = zd z +j, S^ = -d z (2.17) 

and satisfy the usual SL(2) algebra 

[S+,S-] = 2S , [So, 5±] = ±S± . (2.18) 

The quadratic Casimir operator for functions of one variable is just a number: 

S 2 = S 2 - So + , [S 2 ,Si} = 0, S 2 =j(j-1). (2.19) 

The separation of quark and gluon fields in "plus" and "minus" components simultaneously 
serves to separate the terms that transform according to different representations of the 
collinear conformal group. One can show [27, 29] that "minus" projections of the quark 
fields ip—,x- correspond to j = 1/2, "plus" projections of the quark fields ip + ,x + together 

with /_| , /-| gluons transform according to j = 1, whereas /++, /++ have spin j = 3/2. 

Multiplicatively renormalizable quark-antiquark operators of leading twist can be con- 
structed (see e.g. [27, 29]) as highest weights of irreducible components in the direct product 
T J=1 (g> T J=1 corresponding to the total conformal spin j = N + 2 where N is the total 
number of derivatives. E.g. for the left-handed quarks 

o N ( y ) = (-d + ) N My)cT ( ^ + ~^ + ) My) (2.20) 

\D+ + D+J 
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where (x) is the Gegenbauer polynomial, = d^—igA^ is the covariant derivative and 
<9 + = d^n^ = (l/2)X a d a aX a - For simplicity we tacitly assume flavor nonsinglet operators. 

The generators of conformal transformations acting on products of fields at different 
light-cone positions (light-ray operators) are given by the sum of the generators acting on 
the field coordinates, e.g. 

s (n,h) = s (n) + S U 2 ) = z 2 dzi + z 2g z2 + 2jizi + +2j2Z2 

S (Jl,h,h) = z 2q zi + z 2q z2 + z 2q z3 + 2jiZi + +2hz2 + 2hzs ^ (2 21) 

etc. For brevity we will sometimes omit the conformal spin indices if their values and also 
the number of the fields are clear from the context. 

Light-ray operators can be viewed as generating functions for local operators of a given 
twist. Let 

++ (z u z 2 ) = ^+{z 1 )[z 1 ,z 2 ]^ + (z 2 ) , (2.22) 
where ^ 

[z u z 2 ] = PexpfygznJ duA ++ (z^ (2.23) 

is the light-like Wilson line. Here and below we use the following shorthand notations: 

zi 2 = z\ — z 2 , z 21 = uz 2 + uz\ , u = 1 — u . (2.24) 

The conformal operator (2.20) can be written in terms of 0++(zi, z 2 ) as 

N (y) = {d Zl +d Z2 ) N cf f ^ 21 "^ 2 ) ++ ( Zl n + y,z 2 n + y) . (2.25) 

\0 Zl + Z2 J z i= 

Alternatively, one can expand the light-ray operator (2.22) in terms of conformal operators: 
0++( Zl ,z 2 ) = Cdu{uu) N+l O N {z^) 

N J ° 

= E ^ z i2 E Ti f du ^ N+1 ( z 2i) k VNk(0) , (2.26) 



where 



and 



2(27V + 3) 

KN ~ (N + iy. (2 - 27) 



Nk = d« + N . (2.28) 

The operators Ojvfe, k = 0, 1, 2, . . . form what is usually referred to as conformal tower [27] 
and the lowest dimension operator in a tower, Ono = Ctv, is called a conformal operator. 
A formal definition is that a conformal operator is annihilated by the special conformal 
transformations, 50n = [n^K 11 ,On] = (see below). 

Note that the interpretation (or definition) of a light-ray operator as the generating 
function means that one has effectively to deal with polynomials in light-cone coordinates 
Zj. It turns out to be convenient to consider Z{ as complex variables and exploit the 
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equivalence of the SL(2,R) and SU(1, 1) groups. A SU(1, 1) transformation is denned by 
the same formula (2.16), where a,b,c,d are now complex numbers such that d = a and 
c = b: 

. az + b , 
z = ■ (2.29) 

bz + a 

The operator T J , Eq. (2.16), is a unitary operator with respect to the following scalar 
product: 

^ )j = 2 A_± J d 2 z(l-\z\ 2 )^-mz^(z)^ JvjzfoMz), (2.30) 

\z\<l |z|<l 

where the integration goes over the interior of the unit circle (i.e. over the unit disc) in 
the complex plane. In what follows we will often drop the subscript j in the notation 
for the scalar product (. and the measure jT>jz if this cannot yield confusion. The 
SU(1, l)-invariance implies that 

(T^,T>) = (0,^) (2.31) 
which can be checked by explicit calculation using that 

d 2 z'=r, ^T- v>- (2-32) 



(bz + a) 2 (bz + a) 2 
It is easy to show that for simple powers 



<.V> = «W||.l 2 , \\z n \\ 2 = ^^-. (2.33) 

r(n + 2j) 

The SU(1, 1) generators are defined as in Eq. (2.17) and have the following hermiticity 
properties (with respect to the above scalar product): 

Sl = S , (S+? = -S-. (2.34) 

Many of the specific techniques that make this realization useful for our present purposes 
are based on the following representation of the unit operator (reproducing kernel) 



where 



(z) = / VjwlCjfawMw), (2.35) 

|iu|<l 



K ' M = JT=^&fi (2 - 36) 

which will be used in what follows. 

This construction is easily generalized for functions of several variables. The generators 
of conformal transformations are defined as the sums of the generators in Eq. (2.17) with 
spins ji, . . . ,j n which act on the corresponding variables z%, . . . , z n . The scalar product is 
also modified in an appropriate way, e.g. for two variables it takes the form 



( ( t>^)ji,j 2 = J VjiZi J v j2 z 2 <t>{zi, z 2 )ip(zi,z 2 ) ■ (2.37) 

|zi|<l |z 2 |<l 
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One can show that 

n„l|2 r(2j!)r(2j 2 ) r(2n + 2j 1 + 2j 2 -l) 

II uWhh T(n + 2ii)r(n + 2j 2 ) r(n + 2ji + 2j 2 - 1) ' 



With the help of the SU(1, 1) scalar product one can project the conformal operator 
(2.20) from the light-ray operator 0++(zi, z 2 ) as follows [35] 

On = Pn(zi2, ++ (zi,z 2 )) n = Pn J J T>xZ]Vyz 2 zf 2 ++ (z 1 ,z 2 ) , (2.39) 

\z t \<l 

where 

PN = \(N + l)(N + 2)\. (2.40) 

This representation can be proven in the following way. We remind that the operator is 
called conformal if it is annihilated by the generator of special conformal transformations, 
&kOn = [n^K^, On] = 0. The action of the quantum operator n^K^ on the quantum fields 
can be replaced by the differential operator acting on the field coordinates, in particular 

8 k O ++ (z 1 ,z 2 ) = 2(nn) S i +' 1) ++ (z 1 , z 2 ) , (2.41) 

where 

= S+ + S+ = zjd Zl + z 2 2 Z2 + 2*i + 2z 2 . 
Applying the special conformal transformation to Eq. (2.39) one obtains 



SkOn = pn J J T>z 1 Vz 2 z&6kO ++ (z 1 ,z 2 ) ~ J J Vz 1 Vz 2 z? 2 S^ ± >0 ++ (z 1 ,z 2 ) 

\z t \<l \zi\<l 

Vz 1 Vz 2 (s { y ) z? 2 Yo ++ (z 1 ,z 2 ) = 0, (2.42) 



|zi|<l 

as expected. The result is zero because S^ 1,1 ^ = —{d zi + <9 22 )> and thus S_\zi — z 2 ) N = 0. 
It follows that the l.h.s. of Eq. (2.39) is indeed a conformal operator, up to a (convention- 
dependent) coefficient which can be calculated explicitly. A useful formula is 

^11 d 2 Zl d 2 z 2 & - z 2 ) N ^ +i ^ = i N ( Pl + p 2 ) N cf ( . (2.43) 

^ J J \Pl +P2J 

\Zi\<l 

In the same way one finds that 

N k = PN(^ t Nk(zi,z 2 ),0 ++ (z 1 ,z 2 )) u = p N J J V 1 ziViz 2 y t N k 2 (zi,z 2 )0 ++ (zi,z 2 ) , 

N<i 

(2.44) 

where 

^ 2 (z 1 ,z 2 ) = (S ( ^ 1) ) k z? 2 (2.45) 
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can be thought of as a "wave function". The superscript t = 2 serves to remind that we 
are dealing here with twist-two operators. 

The functions ^^(zi, z 2 ) form an orthogonal basis 

<^ fc 2 , ^,|„) = W<W|I^ 2 IP , (2-46) 

so that one can represent the nonlocal operator in the form 

oo oo 

++ (z 1 ,z 2 ) =Y J Y.PN 1 \\^Nkr 2 ^ t Nk^i^2)0 Nk . (2.47) 

N=0 k=0 

The norm 1 1 ^Jv^ 1 1 2 can be calculated recurrently 
One obtains 



t=2 i|2 
Nk-l\\ • 



so that finally 



where 



||**=2||2_ Mr Nu2 n -l „ 1 r(2iV + 4) 

11***11 -ll^lluPM. PNk -k\T{2N + A + k) (2 - 48) 

oo oo 

++ ( Zl ,z 2 ) = ^2^Nk^Nk(^^2)O Nk (0) , (2.49) 

AT=0fc=0 

1 T(N + 2)T(N + 2) 

Explicit expression for the "wave functions" ^^(zi, z 2 ) can be obtained observing that 
«S+ is the generator of special conformal transformations so that expfaS+J corresponds to 
a finite conformal transformation z — > z/(l — az). Therefore 

exp^f ' h) ]z? 2 - ^ 



1- a Zl ) N + 2 n(l -az 2 ) N + 2 » 
T[2N + 2j 1 + 2j 2 } f zg^'-^-i 

/ 1 1 \ 9 /V -1-2-71 ' v * / 



" T[N + 2j x ]T[N + 2j 2 ] J (1 - az t 21 ) 2N + 2 n+^ ' 

where ji = j 2 = 1 for the case at hand. In order to find ^^(zi, 2:2) one has to differentiate 
this expression k times in a and set a = 0. This gives: 

r S (n>hhk z N = N T[2N + 2j 1 + 2j 2 + k] f 1 N+2jl - HN+2h -! , t ) k (2 _ 52) 

{+ ' 12 12 r[iv + 2ji]r[iv + 2j 2 ] J 1 2i; 1 ; 

Using this representation the expansion (2.49) can easily be recast in the form (2.26). 
What has to be stressed here is that the functions which project the local operators from 
the nonlocal one are the same functions that appear in the expansion of the nonlocal 
operator over the conformal operators and its descendants (conformal tower). This is a 
general property which holds for arbitrary operators [29, 36]. 
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2.3 Conformal non-quasipartonic operators 

All existing (gauge-invariant) composite operators in QCD can be divided in two classes: 
quasipartonic [28] and non-quasipartonic. By definition, operators built of "plus" compo- 
nents of the fields and arbitrary amount of "plus" covariant derivatives are called quasi- 
partonic. The light-ray operator defined in Eq. (2.22) is the generating function for local 
quasipartonic operators of the leading twist. The complete operator basis for the leading 
twist-two and also subleading twist-three operators can always be chosen to contain quasi- 
partonic operators only. The subset of quasipartonic operators of a given twist is special 
in that it is closed under renormalization [28] at one-loop. Moreover, since building blocks 
of quasipartonic operators transform as the primary fields under SL(2, M) transformations, 
the corresponding RG equations are explicitly SL(2, M)— invariant [28, 29]. Examples of 
twist-four quasipartonic operators are V>+'0+V'+V'+7 ip + f ++ f ++ ip + , etc. 

All other operators are naturally called non-quasipartonic, e.g. they can include "mi- 
nus" fields components and "minus" or transverse derivatives. Non-quasipartonic operators 
cannot be avoided starting with twist-four. The complete operator basis for twist-four op- 
erators, including non-quasipartonic ones, was constructed in [29]. 

The main problem in dealing with non-quasipartonic operators is due to transverse 
derivatives which, generally, do not have "good" conformal properties. Note that in the 

spinor formalism we distinguish two transverse derivatives: and D (_. One can 

show [29] that e.g. ip + transforms as a primary field with conformal spin j = 3/2, 

whereas D \ r ip + does not transform homogeneously, it is a "bad" object as far as the 

SL{2) symmetry is concerned. The procedure suggested in Ref. [29] is to keep the "good" 
derivatives and get rid of the "bad" ones using equations of motion D aa ip a = etc. The 
possibility to make this separation is the crucial advantage of the spinor formalism in the 
present context. 

In this work we will be interested in operators with quantum numbers such that they 
can contribute to the time-ordered product of two electromagnetic currents. With the fur- 
ther restriction to flavor- nonsinglet contributions, we are left with three (non-quasipartonic) 
operators of twist-four that involve the self-dual gluon field /, 




and three operators involving /: 




The rationale for including the factors (/xA) and (A//) in the definition of the operators is 
that, in this form, they contain equal amount of auxiliary spinors with and without a "bar" 



(i.e. equal amount of A and A and separately fj, and fj) and therefore can be rewritten in 
terms of auxiliary vectors n, h, e.g. 

Qi(n, n; z\, z 2 , z 3 ) = $(zin)nf(z 2 n)nn-if)(z 3 n) . (2.55) 

Note that 

(/xA)(A/2) = 2(n-n). 
Operators in usual vector notation can easily be converted to this basis, e.g. 

(nn)q L (zi)[F +fl (z 2 ) + iF +hl (z 2 )]^q L (z 3 ) = Q 2 {zi,z 2 ,z 3 ) - Qi(zi,z 2 ,z 3 ) , 
(nh)q L (z 1 )[F + ^ t (z 2 ) - iF +fl (z 2 )]^q L (z 3 ) = Q 2 (zi,z 2 ,z 3 ) - Qi(zi,z 2 ,z 3 ) , (2.56) 

where q L = |(1 - 75)9. 

The operators with and without a "bar" are related by hermitian conjugation: 

Qi(z 1 ,z 2 ,z 3 ) = (Qi(z 3 ,z 2 ,z 1 ))' f (2.57) 

and transform as direct products of primary fields under SL(2, R) transformations, as indi- 
cated. We stress that the operators that include a transverse derivative, Q 3 and Q 3 , must 
be included in the basis in order to ensure simple realization of the conformal symmetry, 
although they may be dispensed off at the later stage (see below). 
The following combinations 

Qk(zi,z 2 ,z 3 ) = Qk(z!,z 2 ,z 3 ) ± Qk(z 3 ,z 2 ,zi) (2.58) 

transform in a different way under the combined charge conjugation and parity trans- 
formations and therefore cannot mix under renormalization. By a direct calculation one 
finds 

4- CP 4- 

Q k (n,h;zi,z 2 ,z 3 ) — > ±Q k (n, n; z\, z 2 , z 3 ) , (2.59) 

where we have displayed explicitly the dependence on the auxiliary vectors n and n. Note 
that the auxiliary vectors are interchanged nf>nby tbe CP transformation. For compar- 
ison, the leading-twist light-ray operator (2.22) transforms as 

CP 

++ (n, z 1 ,z 2 ) — > -0 ++ (h, z 2 , zx) , (2.60) 

which implies that for conformal operators (2.39) 

N (n)^{-l) N+1 N {fi). (2.61) 

Note that interchanging the auxiliary vectors does not have any physical significance as 
they only serve to simplify the leading-twist projection for local operators. 

To save space, in the following discussion we will sometimes use vector notation: 

/ Qi(zi,z 2 ,z 3 ) s 
Q 2 (zi,z 2 ,z 3 ) I (2.62) 



Q{z 1 ,z 2 ,z 3 ) 
and similar for Qi. 



\Q3(zi,z 2 ,z 3 ) 
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2.4 Renormalization and hermiticity 

Let P(zi, z 2 , Zs, Z4) be the complete set of quasipartonic twist-four operators of the type 
, ijj + ip + ip + ip + , -ip + f ++ f ++ ip + , etc. Explicit expressions will not be needed in what follows. 
The complete RG equation (to the one-loop accuracy) has a block-triangular form 



9 0< N 9 




a s 
2^ 










M 


i 






Q 


V 





Upp) 





(2.63) 



so that the quasipartonic operators have autonomous evolution, whereas the non-quasi- 
partonic ones, Q and Q , do not mix with each other but can mix with the quasipartonic 
operators. Due to this structure, the RG equation (2.63) actually decouples in two equa- 
tions which contain Q and Q operators: 



9 n / \ 9 

9 Q { \ 9 

^ + /?(as) 9a- 





QQ ^QP 
Upp 

QQ ^QP 
Mpp 




(2.64) 



The operators Qi and Qi are hermitian conjugate to one another, cf. (2.57), so that the 
two RG equations in (2.64) are equivalent and it is sufficient to consider one of them. 
This implies in particular that the spectrum of one-loop anomalous dimensions of non- 
quasipartonic operators is double-degenerate: there are two operators with different CP- 
parity for each anomalous dimension. 

The "Hamiltonians" EI are known explicitly [30] and can be written in terms of SL(2) 
invariant integral operators (kernels). "Diagonal" Hamiltonians Hqq, M.qq and Mpp have 
a two-particle structure, i.e. each contributing kernel involves light-cone coordinates of two 
field operators only, whereas the mixing ones, Hqp and Hq p , involve 2—7-3 transitions. 

The light-ray operators , Q can be expanded in multiplicatively renormalizable 

local operators, cf. Eq. (2.47). Thanks to the block-triangular structure of the mixing ma- 
trix, expansion of ~P only involves quasipartonic local operators, which we denote by Vnp, 
whereas the expansion of Q (Q) involves both Vnp and multiplicatively renormalizable 
non-quasipartonic operators, which we denote by Qn p (Q.Np)' 



3 



Q 



^ aNpjzi, z 2 , z 3 , z^)V Np , 
Np 



^ bN P (zi,z 2 , z 3 )Q Np + ^2 CNp(zi,Z 2 ,Z 3 )V Np , 
Np Np 

y~] b~Np{zi,z 2 , z 3 )Q Np + ^ CNp(zi, Z2, Z 3 )V N p , 
Np Np 



(2.65) 



where N is the total amount of derivatives (which specifies the operator dimension) and p 
is a multi-index that enumerates different operators. 
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Substituting this expansion in the RG equation (2.63) one finds that the coefficients of 
non-quasipartonic local operators 6jvp(^i, Zz, £3) (pNp{zi, %2, Z3)) can be found as solutions 
of the integral equation that only involves Mqq 



(2.66) 



lb Np (z 1 ,Z 2 ,Z 3 ) = jNpbNp{zi,Z2,Z 3 ) , 
t b Np (zi,z 2 ,z 3 ) = jNpbNp(z 1 ,Z 2 ,Z 3 ) , 



where 7at p are the corresponding anomalous dimensions. The two equations in Eq. (2.66) 
are related by hermitian conjugation, cf. (2.66), so that l)Np(zi, z 2 , z 3 ) = ±bN p (z 3 , z 2 , z\), 
depending on the parity of the operator Qn p under the CP-conjugation. 

In the context of this work we are interested in a particular set of contributions that 
are related to the divergence of the leading twist conformal operator On, which we denote, 
schematically, (80)n, and its descendants obtained by adding arbitrary amount of "plus" 
derivatives, d\{dO)N. These operators can be identified as the non-quasipartonic twist- 
four operators with the anomalous dimensions coinciding with the anomalous dimensions of 
leading- twist operators. Hence the corresponding coefficient functions b^kizi, z 2 , z 3 ) corre- 
spond to the subset of solutions of the integral equation in Eq. (2.65) with the eigenvalues 
7at that are usual flavor-nonsinglet leading-twist anomalous dimensions: 



JV+l 



IN = C F 1 



1 



(N + l)(N + 2) 



+ 4V- 



m=2 



in 



2C F 



iP(N + 3) + i/)(N + 1) - V(3) - -0(1) 



(2.67) 

Note that the quasipartonic operators completely decouple and their contributions can be 
ignored at any stage of the calculation. 

The operator Mqq is a 3 x 3 matrix with the entries being two-particle integral oper- 
ators: 



Hn H\2 Hi 3 
H 2 i H 22 H 23 
H31 H 32 H 33 



(1 ^ Al 

(l.l.l) 



(2.68) 



where the superscripts on the three components of the "wave function" indicate the corre- 
sponding SL(2) representations. Explicit expressions [30] are given in Appendix A. 

Unfortunately, Eq. (2.66) appears to be too complicated to be solved directly. The 
way out is that the operator Mqq turns out to be self-adjoint (hermitian) with respect to 
the following scalar product: 



((¥,¥; 



2(*1, *l)lll + ($2, ^2>i3 1 + -($3, 

22 2 



3/33 1 , 
2 2 1 



(2.69) 



which can be verified by the explicit calculation, cf. Appendix A. 

Thanks to hermiticity it proves to be sufficient to solve, instead of Eq. (2.66), the 
(much simpler) inverse problem: find the expansion of the divergence of the leading-twist 
conformal operator {80)n in terms of quark-gluon non-quasipartonic operators. Examples 



-17- 



of this expansion for N = 1,2 are given in the Introduction, cf. Eqs. (1.2), (1.3). We have 
been able to obtain this expansion for arbitrary N and bring the answer to the form 



where the functions 



WPn 



^n(zi,z 2 ,z 3 ) 



A,3)>-«*AT,5>> 



+ 



(2.70) 



(-l) N ~$ N (z 3 ,Z 2 ,Z 1 ) (2-71) 

are known explicitly and the ellipses stand for (irrelevant) contributions of quasipar- 
tonic operators. The corresponding calculation is presented in Sect. 3. The symmetry 
relation (2.71) follows from properties of the operators under CP-transformations, see 
Eqs. (2.59) and (2.61). 

Note that the two contributions in the square bracket in Eq. (2.70) are renormalized 
multiplicatively with the same anomalous dimension 77V (2.67). Thus we can construct 
another multiplicatively renormalizable twist-four operator 



(juA)(A//)Tjv 



(N + l) 



+ 



(2.72) 



which has the same leading twist anomalous dimension and is not reduced to total deriva- 
tives of the conformal operators. Hence it contributes, e.g., to the total cross section of 
deep-inelastic scattering. An example will be given below. 

The same expansion holds for the full conformal tower obtained by adding "plus" 
derivatives to (dO) N , cf. Eq. (2.44): 



(fiX)(xp)di(do) N 



wpn 
(N + l) 



Us + )*# N ,Q))-(({s + y* N ,Q)) 



+ ... 



(2.73) 



where the operators S^. 1 '^ 2 = z\d zx + z\ d Z2 + z\d Z3 + 2j\Z\ + 2j 2 z 2 + 2j 3 z 3 acting on the 
three components of ^ n have to be taken in the corresponding representations, i.e. 



5+* 



I 5(1,1,1)^(1,1,1) ^ 

g\ ' 2 ' 2 > 



(3 3 ^ , 3 3 jn 
' 9 - \J/ 2 ' 2 



(2.74) 



The functions ^ j\rfc = (S + ) are mutually orthogonal, and also orthogonal to the 
coefficient functions of the other existing multiplicatively renormalizable operators. It is 
easy to show that 

P^ll^vll 2 (2-75) 



|*JVfc| 



with the same coefficient as for the two-particle leading-twist operators, Eq. (2.48). It 
follows that the light-ray non-quasipartonic operators £^ can be written as 



ifj 



zi,z 2 ,z 3 ) 



OO OO , ,j ^ s 9 

W E E IUT, 12 ^Nk(z!,Z 2 , Z 3 ) 8%(dO) N + ... 

pn\\^n\\ z 



igQ(z!,z 2 ,z 3 ) 



N=l k=0 

00 00 

w E E 

N=l k=0 



p Nk (N + lf 

A/v||*!v|| 2 



(2.76) 

l) N+1 ^ Nk (z 3 ,z 2 ,z 1 )d k 2(dO) N + ... 



-18- 



which is the desired result. The ellipses stand for the contributions of all other exist- 
ing twist-four operators that cannot be reduced to the "plus" derivatives of (80)n- As 
mentioned above, one of such "genuine" twist-four contributions has the same anomalous 
dimension as in the leading twist. 

2.5 Reduction of the operator basis 

The description of the contribution of "kinematic" operators d+(dO)N in terms of the 
three-component "wave function" \f = {^i, Vf^j ^3} as in Eq. (2.77) is in fact redundant. 

As discussed in [29], the expansion of light-ray operators in terms of local multiplica- 
tively renormalizable operators (2.65) runs over all operators of collinear twist-4. Thus 
together with the operators of geometric twist-4, this expansion contains the operators 
which are descendants of the geometric twist-3 operators. 2 The latter ones are not inter- 
esting in the present context since they have autonomous scale dependence and decouple 
in the similar way as the quasipartonic geometric twist-4 operators. 

It can be shown (see below) that the functions n p which correspond to the operators 
of geometric twist-4 satisfy the following conditions: 

— <9 Z2 4i*i (21,^2, z 3 ) + d Z3 Z3i^ 2 {zi, Z2,zs) = 0, (2.77a) 

—d Z2 Z^i(z 1 ,Z2,Z 3 ) + ^2(z 1 ,Z2,Z3)+Zi3^3(zi,Z2,Z 3 ) = 0. (2.77b) 
Z23 

These relations are SL (2)-invariant, i.e. if Eqs. (2.77) hold for the components of n p , 
they hold for the functions 5±,o * n p as well, where it is assumed that the generators are 
taken in the corresponding representation, cf. Eq. (2.74). 

The relations in Eq. (2.77) hold for the coefficient functions of an arbitrary operator of 
geometric twist-4, and in particular they hold for the coefficient functions of the divergence 
of the leading twist operators, Eq. (2.77). In principle, one can resolve these equations and 
express two of the functions in terms of the third one, e.g. and ^3 in terms of 

rzi 

z 21 ^/ 1 (zi,z 2 ,z 3 ) = d Z3 z l3 dw (zi - w)^ 2 (zi,w,z 3 ) , (2.78) 

J Z2 

however, using the "three-component" formulation may be more convenient because the 
scalar product (2.69) has a simple form. 

In order to derive Eq. (2.77), consider the quasipartonic twist-3 light-ray operator 
Qs(zi, z 2 , z 3 ) = ■ijj + (zi)f ++ (z 3 )ijj + (z 3 ). Similar to the above, Q 3 can be expanded in con- 
tributions of multiplicatively renormalizable local operators: 

Q 3 (z u z2,z 3 ) = J2^~ 3 (zi,z 2 ,z 3 )Q^- 3 . (2.79) 

Nq 

2 We remind that geometric twist of an operator is denned as "dimension minus spin" whereas collinear 
twist is given by the difference in dimension and spin projection on the light cone. The classical example of 
the contribution of geometric twist-2 operators to the collinear twist-3 observable is the Wandzura-Wilczek 
contribution [3] to the structure function g2(x,Q 2 ) in the polarized deep-inelastic scattering. 



- 19 - 



The trick is to calculate the derivatives fj, a (d/d\ a )Q 3 (z) and [iP^\, Q 3 {z)] in two different 
ways: On the one hand, making use of the expansion in Eq. (2.79) one obtains 



-Nq 



fid x Q 3 ( Zl ,z 2 ,z 3 ) = ^^ 3 (z 1 ,z 2 ,z 3 )^AQf" 

Nq 

\iP u - x Qs(z 1 ,z 2 ,z 3 )\ = y^%-\z 1 ,z 2 ,z 3 )\iP fl - x ,Q%- 3 } 



Nq 



(2.80) 



The sums on the r.h.s. of (2.80) do not contain genuine geometric twist-4 operators, but 
rather the descendants (derivatives) of the twist-3 operators only. On the other hand, one 
can express [id\Q 3 (z) and [iP^x, Qs{z)] directly in terms of the (^-operators: 

(\p)[idxQs(zi,Z2, z 3 ) = z 2 1 d Z2 z 2 Qi(z 1 ,z 2 , z 3 ) + d Z3 z 3 Q 2 ( Zl ,z 2 ,z 3 ) + z x Q 3 (zu z 2 , z 3 ) , 
(^M^Va^^i,^,^)] = d Z2 Qi{zi,z 2 , Z3 ) + d zs Q 2 { Z i, Z2 ,z 3 ) + Q 3 ( Z i, Z2 ,z 3 ). (2.81) 

Since the expansion of Qk over the local operators (2.65) contains all operators of collinear 
twist-4, it means that the geometric twist-4 operators must drop out in the particular 
combinations appearing in (2.81) and only the descendants of geometric twist-3 operators 
survive. If follows that the coefficient functions of the geometrical twist-4 operators in this 
expansion must satisfy the equations 



{z 2 d Z2 +2)*i(2i, z 2 ,z 3 ) + (z 3 d Zs + l)^ 2 (zi,z 2 ,z 3 ) + zi^ 3 (zi, z 2 , z 3 ) = 0. 

^2*1(^1,22,^3) + ^3*2(^1,^2,23) + $ 3 (z 1 ,z 2 ,z 3 ) = 0, 
which are equivalent to Eqs. (2.77). 



(2.82) 



3 Divergence of a conformal operator 

Let 0^ Ll "' tJ ' N be a multiplicatively renormalizable leading-twist operator with N covariant 
derivatives, symmetric and traceless over all open indices. The conformal operator On 
defined in Eq. (2.20) is obtained by the projection 

O n = nV . . . n^{0 N ) mi ... m , (3.1) 

where n is an auxiliary light-like vector. We define a divergence of the conformal operator 



as 



(00) jv = ...n» N d»(0 N ) 



n 



/'i 



n 



I'N 



i*" 1 , (On) 



where P. is the usual four-momentum operator 



Pn\. 



d 



$(x) . 



Taking into account that = i(Ao"^A) the same definition can be rewritten as 



(dO) 



1 



N 



N + l 



1 



(3.2) 



(3.3) 



(3.4) 



(iV + 1) 2 

where the argument of Oat indicates that it has to be considered as a function of or 
A Q ,A°, respectively. 
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3.1 Expression in terms of non-quasipartonic operators 

Our first task is to rewrite this expression in terms of quark-gluon operators using QCD 
equations of motion (EOM). To this end it is convenient to use the integral representa- 
tion (2.39) for the conformal operator On in terms of the light-ray operator 

++ (z 1 ,z 2 ) = 'ip + (z 1 )[z 1 , z 2 }ip + (z 2 ) = 'tp(z 1 n)fi[z 1 n, z 2 n]'ip(z2n) . (3.5) 

A subtlety is that this operator is defined on the light-cone n 2 = and in order to take 
derivatives with respect to one has to extend this definition to n 2 ^ 0. This can be 
done by applying the trace subtraction operator (see the next section) , but a more elegant 
way out is to use the last (spinor) representation in Eq. (3.4). Indeed, taking derivatives 
with respect to the auxiliary spinors A, A one stays on the surface n 2 = 0. 
Let 



(dO) ++ ( Zl ,z 2 ) 



(3.6) 



so that 
(dO) N 



PN 



z? 2 ,(dO) ++ ( Zl ,z 2 )) 



(N + l) 



PN 



(N + 



!) 2 II\ Zi \<i 



Vz{Dz 2 4 (0O) ++ 

(3.7) 

(3.8) 



Taking into account that 





d 
3A« 



one obtains 

(dO) ++ (z u z 2 ) 



iP„,((nd) + l)0 ++ (z 1 ,z 2 ) 



(iPn),0 ++ (zi,z 2 ] 



(3.9) 



where = Making use of the identities 



{(nd) + l)0++(zi, z 2 ) = ( Zl 8 Zl + z 2 8 Z2 + 2)0++(zi, z 2 ) 
(iPn),0 ++ (z 1 ,z 2 ) = (d zi + d Z2 )0 ++ (z 1 ,z 2 ) 



this becomes 
(dO) ++ (z u z 2 ) 



zid Zl + z 2 d Z2 + 2) 



(3.10) 



iP^d^o ++ ( Zl ,z 2 )\ - -(d Zl +d Z2 )d 2 o ++ ( Zl ,z 2 ) 



(3.11) 

Note that the definition in Eq. (3.6) involves the light-ray operator ++ (zi, z 2 ) at strictly 
light-like separations n 2 = and the calculation leading to Eq. (3.11) is pure algebra. 
Hence the expression in (3.11) effectively does not contain derivatives in the directions 
"orthogonal" to the light-ray, so that one can treat n^ 1 as a generic four- vector and put 



n 



at the end of calculation. 
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The derivatives of the light-ray operator in Eq. (3.11) include contributions from the 
Wilson lines [37, 38]: 



iPfj,, [zin,z 2 n] = igA^(zi)[zin, z 2 n\ - ig\z\n, z 2 n]A^(z 2 ) 



(3.12a) 



dn, 



+ igz 12 / du[zin,Z2 1 n]n u F^ u (z 21 )[z 2l n,z 2 n] 
Jo 

[zin,z 2 n\ = igz 1 A fJ ,(z 1 )[z 1 n, z 2 n] - igz 2 [zin, z^A^zi) 



+ igz\ 2 / du z 21 [zin, z\\v\n v F^ v {z\^\z^ x n, z 2 n] 



(3.12b) 



In particular, taking into account Eq. (3.12b) one obtains 

d^O ++ (zi, z 2 ) = 4>{zin)d^il){z 2 n) + z 1 [D^}(zin)nip(z 2 n) + z 2 ${zin)fi[D ll ip](z2n) 



igz\ 2 \ du Z2i^{z\n)nn v F^ y {z2i)^{z 2 n) 



(3.13) 



where for brevity we do not show the Wilson lines in the operators on the r.h.s. The 
covariant derivatives are defined as D^ip = (d^ — igA^t 11 )^, = (c^ + igA a ^{t a ) T ) r il). 

Calculating the commutator of (3.13) with iP^ (or taking the second derivative d^) 
one gets contributions of the following type: 



D^nD^ip, D 2 ipnip, tpnD 2 ip, igipFvp , ipFFip, ipD^F^tf; . 



(3.14) 



The last term ~ i^D^F^ip can be reduced to the four-quark quasipartonic operator using 
EOM. Also the terms ~ ipFFip are quasipartonic. As explained above, such contributions 
are irrelevant for our discussion and will be omitted. 

The contributions of the first type, ~ D^ipnD^ip, are, on the other hand, the only 
ones that do not vanish in the free theory. One should expect that such terms do not 
contribute to a divergence of the conformal operator by virtue of the FGPG theorem [20]). 
It is instructive to see how this result is recovered in our formalism. 

The contribution to {dO) ++ {z\, z 2 ) (3.11) due to the terms ~ D^nD^tp takes the 
form 



1 



(dO)f;( Zl ,z 2 ) = \[ Zl d zl + z 2 d Z2 + 2)(zx+z 2 ) - -(d Zl +d Z2 )2z lZ2 \ [D^\{ Zl n)n[D^](z 2 n) 



{z\d Zl + z 2 2 d Z2 + 2 Zl + 2z 2 )[D^]( Zl n)n[D^](z 2 n) 
S+[D^]( Zl n)n[D^}(z 2 n), 



(3.15) 



so it is not zero, but proportional to the two-particle "step-up" operator (2.17). It follows 
that 

(dO)% ee ~ (z^St 2 [D^]( Zl )n[D^](z 2 )) = -<Sr 2 4, [D^]( Zl )n[D^](z 2 )) = 0, (3.16) 



as expected. 
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The remaining terms ~ D 2 ipnip, ~ ipnD 2 ip and ~ ig^Ftp give rise to the contributions 
of the non-quasipartonic operators of interest (2.53), (2.54). The corresponding calculation 
is detailed in Appendix B. The answer has the form 



(n\){\n){dO) ++ {zi,z 2 ) = ig A(zi,z 2 ) - A(zi, z 2 ) 



+ ... , 



(3.17) 



where 



A(z 1 ,z 2 ) = d Z2 zf 2 <Q 1 (z 1 ,z 1 ,z 2 ) + / duu Q 2 (z 1 , z 21 , z 2 ) + z 12 Qz{z\, z 21 ,z 2 ) 



+ d Zl d Z2 zf 2 du - Qi {zi,z 21 , z 2 ) +u Q 2 (z 1 , z 21 ,z 2 ] 



(3.18) 



A(zi,z 2 ) = d Z2 z 21 \Q 1 (z 1 ,z 2 ,z 2 ) + I duu Q 2 (z 1 , z" 2 , z 2 ) + z 2 iQz{z\, z± 2 , z 2 ) 



+ d zl d Z2 z 21 du - Qi(zi,zi 2 ,z 2 ) +uQ 2 (z 1 ,Zi 2 ,z 2 ) 



(3.19) 



The ellipses stand for EOM, contributions of quasipartonic operators and terms propor- 
tional to S± 2 which do not contribute to the projection (3.7) that defines the conformal 
operator: 



(ji\)(\p)(dO) N 



3.2 Coefficient functions 



{z 12 ,A(zi,z 2 ))u - (z 12 ,A(zi,z 2 )) n 

Vz{Dz 2 z^ 2 [A(z 1 ,z 2 )-A(z 1 ,z 2 )] . (3.20) 

«fe|<i 



Our answer for the divergence of the conformal operator (3.20) involves the sum of terms 
of the type, schematically 



(aO) N ~{z^ i Q i ](z 1 ,z 2 )) 11 , 



(3.21) 



where Bj are some integral operators, and the subscript (, )n indicates that the scalar 
product is calculated for conformal spins j\ = j 2 = 1. The idea is to rewrite this answer 
as a sum of terms 

5 z 2i z 3))jij2j3 ' (3.22) 

where ji,j 2 ,j3 are the conformal spins of the Qi, so that the functions ^f{z\, z 2 , z%) can 
be identified with the coefficient functions of the Q-operators, cf. Eq. (2.70). 

For illustration let us consider contribution of the first two terms in (3.18) which we 
write as 

A = d Z2 zl 2 [0i (z u z 2 ) + <fa(zi,zz) + ..-] + ... , (3.23) 

where 



01 {zi, z 2 ) = Qi(z 1 ,z 1 ,z 2 ) , 
(^2(^1,22) = / duuQ 2 (zi,z 21 ,z 2 ) 



(3.24) 
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As the first step, notice that any operator Qi(z±, Z2, 23) can be written as 

3 

Qi(z 1 ,z 2 ,z 3 ) = (Y[fcj k (zk,w k ) ) Qi(wi,W2,w 3 )). i . 2j3 , (3.25) 

k=l 

where lCj k are the reproducing kernels defined in Eq. (2.35) and the scalar product corre- 
sponds to the spins j\ , j'2 , j 3 . Thus 

<Pl(zi,z 2 ) = (fci(zi,wi)lCi(z 1 ,W2)lC 1 (z2,W3)Qi(wi,'W2,w 3 )) 111 , (3.26) 

<£2(2l,22) = / duu(K:i(zi,Wi))C 3 /2(z%l,W2)lCi/2(z2,W3)Qi(wi,W2,W 3 )) 31 

Jo 22 
= -(/Ci(zi,U)i)/Ci(zi, W 2 )JC 1 / 2 (z2,W2)]Ci/2(z2,W 3 )Qi(wi,W2, Wz)) 1 3i (3.27) 

where we used that 

-1 f 1 1 f 1 u 

duuKy 2 {z 2 n ^2) = / duu- — — 7^3= / du TTi ^ i -n i ^3 

Jo {1 - {z 2 u + ziu)w 2 y Jo {u(l - ZiW 2 ) + u(l + Z 2 W 2 )) 6 

= - jz — 2-rj z-t = -Ki{zuiD 2 )1C l / 2 {z2,W2) . (3.28) 

2 (1 - Z1W2Y (1 - £2^2) 2 ' 

Making the conformal transformation z —¥ z', w — > w 1 and taking into account that the 
reproducing kernels become 

Kj(z',w') = (bz + a) 2j K.j{z, w)(bw + a) 2j , (3.29) 

it is straightforward to check that both (fi(zi,z 2 ) an d f2{zi,z 2 ) transform according to 
the representation T 31 ® T 32=1 . Using this result, it is easy to show that the both 
contributions in Eq. (3.23), d2z\ 2 l Pi{zi, Z2) and 82 zf 2^2(21, Z2), transform according to the 
representation T 31 ® T 32=1 , as they should. 

The transformation properties follow immediately from the following statements which 
can be checked by a direct calculation: 

• if a function ip(zi, Z2) transforms according to the representation T n ® T 32 , then the 
function tp(zi, z 2 ) = z\ 2 ip(zi, z 2 ) transforms according to the representation T 3l ~ l l 2 ® 
7 1 J' 2-1 / 2 ) i. e . multiplication by Z\ 2 = z i — ^2 intertwines the representations 

z 12 T jl ® T h = T jl ~ 1/2 ® T J2 ~ 1/2 0i 2 (3.30) 

• if a function (p(z) transforms according to the representation T 3=0 then its derivative 
d z (f(z) transforms according to the representation T 3 , i.e. 

d z T j=0 = T j=1 d z . (3.31) 

We have verified that all contributions to the functions A(zi, z 2 ) (3.18) and A(z\, Z2) (3.19) 
transform separately according to the same representation, T 31 <g> T 32=1 , which provides 
a strong check of the calculation. 
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The contribution of ^1,2(^1 , ^2) to (80)n has the form 

(dO) N ~ (z? 2 ,Mip( Zl ,Z2))ii, (3.32) 

where the operator B = d 2 z\ 2 intertwines the representations IT 31 ® yj2=i = yji =1 ^ 
T^ 2=1 B, i.e. 

B Sgi =2j ' a=1) = 5^ =lj2=1) B . (3.33) 

As the second step, we want to rewrite (3.32) in the form 

(zi 2 ,Mipi(z 1 ,z 2 ))j 1 =j 2 =i = (M^z^ 2 ,ipi(z 1 ,z 2 )}j 1 =2,j 2 =i- (3.34) 

It is easy to see that B^z^ = b^z^ 1 where 6/v is a numerical coefficient which can be 
fixed as follows: 

(z^Mz^n = -(N+ l)||4llii = (Bt4,4- 1 } 21 = bxWzg-Xi ■ (3-35) 
Taking into account (2.38) one gets 

b N = ~N{N + 2)(N + S), (3.36) 
b 

so that 

{z±2, d 2 zf 2 Pi(zi, z 2 ))u = b N (z^ 2 ~ 1 ,Lp 1 (z 1 ,z 2 )) 2 i = ($l a {Wi,W2,W 3 ),Qi(wi,W 2 ,W 3 )) 1 il , 

(3.37) 

where 3 



^i a (wi,w 2 ,w 3 ) = b N 



// V 2 ziT>iz 2 z± 2 1 K. 1 (w 1 ,zx)}Ci(w2, zi)JCi(w 3 , z 2 ) . (3.38) 

J J\z k \<l 

Combining the two reproducing kernels with the help of Feynman parametrization 

/Ci(t«i, zi)Ki{w2, Z\) = 6 / daaa)C 2 (wf 2 ,zi) (3.39) 

J 

and using that 

f [ V 2 z 1 V 1 z 2 ^-^M. 22) = « 2 - ^s)^ 1 (3.40) 
JJ| 2fe |<i 

one obtains 

Via(w 1 ,w 2 ,w 3 ) = —N(N + 2)(N + 3) / da aa « 2 - ws)^" 1 • (3.41) 

J 

Similarly 

(z^,(9 2 ^i2^2(^l,^2))ll = bjv(^~ 1 ,¥'l('Zl,22))21 = (*2o(l0l,W2>W3),<92(«'l>W2,W3))lll , 

(3.42) 



3 The subscript la serves to remind that this is only one of the two contributions to the coefficient 
function of the Qi operator; the second contribution comes from the term in Qi in the second line in (D.l). 
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where 



V 2a (w 1 ,w 2 ,w 3 ) = --N(N + 2)(N + 3) daaa d/3(wf 2 - w^) 



2 .10 Jo 



J \N-1 



(3.43) 



All other contributions to (3.18), (3.19) can be managed in the same manner. Collecting 
all terms we obtain the final result: 



(zi 2 , A(zi,z 2 ))n = ((^ n(w 1 ,W2,w 3 ),'^(w 1 ,w 2 ,w 3 ))) 



(3.44) 



where the invariant scalar product ((. . .)) is defined in Eq. (2.69) and the coefficient func- 
tions ^ n = |^jv^' ^A^' are gi ven by the following expressions: 



^n\ w ) = 0>N T> 3 / 2 Z 1 T> 3 / 2 Z 2 Z12 1 ^1 (™1 , ^l)^l/2 (™2, Zl)K-l/2 (w 2 , Z 2 )K\ (w 3 , Z 2 ) 



+ 1 1 ' D 2ZiT> 1 z 2 z 12 K, 1 {w 1 ,zi)K, 1 {w 2 ,z 1 )K, 1 (w 3 ,z 2 ) 



(3.45) 



4a N 



dB B (wf 2 — w^ 2 ) N 1 — / 



N-l 



Wn\w) = -a N 1 1 V 3 / 2 ziV 3/2 z 2 z^ 2 1 Ki{wi, z 1 )K l/2 {w 2 , zi)Ki{w 2 , z 2 )Ki /2 {w 3 , z 2 ) 
+ ^b N II V 2 ziV 1 z 2 z^ 2 ~ l K 1 {wi,z 1 )K 1 (w 2 ,z 1 )K, l / 2 {w 2 , z 2 )K, 1/2 {w 3 , z 2 ) 



-Aa N J da a J dB (b + ^ _^ aj (wf 2 - w§ 2 ) 



13 sN-1 



(3.46) 



Here 



-CAT 



J J ^5/2^1^3/2^2 Z^ 2 2 K, 3/2 (wi,Zi)ICi(w 2 ,Z 1 )ICi/ 2 (w 2 ,Z 2 )ICi(w 3 ,Z 2 ) 



24cn / da a 2 a / dB 8 (w 



J 12 ~ w 3 2 



\N-2 



1 



a N 
c N 



(N + 3)(N + 2)(N + 1)N, 



6 



1 



(N + 3)(N + 2)N , 
(iV + 4)! 



(3.47) 



(3.48) 



48 (N + l)(N-2)\ ' 
We have checked that these expressions satisfy the relations (2.77). 

3.3 Calculation of the norm ||\I/jv|| 2 

To pick up the contribution of the operator (80)n to the light-ray operators Qk (2.77) we 
need to know the norm of the "wave function" Vw: 



l^ll 2 = 2\\^ 



N- 

W||2 + ||M/( 2 )||2 + 
N Mill 1 1 JV 1 1 1 3 i "r 



1 lltl/ (3) ll 2 



(3.49) 
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Consider the last term in (3.49) as an example. To evaluate 



\ N ' N / § § 1 



(3.50) 



we use the parametric representation for "ket" function (second line in (3.47)), and the 
reproducing kernel representation for the "bra" function (first line in (3.47)), respectively. 
Rewriting the product of reproducing kernels as 



JC 1 (w 2 ,Zi)IC 1 /2(w2,Z2) = 2 / (IUU1C 3 /2(W2,Z% 1 ) 

Jo 



(3.51) 



and using Eq. (2.35) one can easily take the integrals over wi,W2,W3 to get 



^^lllai = 48c ^ll 2 i2 -2 |ll3 / daa' z a I d/3 (3 I du u (a + u{a - 0)) 



o 



2 1 1 „iV — 2 1 1 2 

2 2 

1 1 1 2[i/>(N + 2) - ^(2)] 



o 



7V-2 







& a N 1 1^12 1 111 



3 2(iV + l)(iV + 2) 
The other two contributions are calculated similarly: 



N(N + 3) 



(3.52) 



|^H 2 = 4^11411?! 



|^ ) || 2 =4a JV ||4ll 1 i 



(iV- + l)(i>T + 2) _ _5_ 

N(N + 3) L ^ v ' ^ v n 12 
2[i/>(N + 2) - ^(2)] 



(N + 3)N 



+ 1 + 



1 



2(iV + l)(iV + 2) 



Collecting all terms we obtain 



l^vll 2 = ^||4lln (N + 2) 2 (iV + l) 2 ty(N + 3) + + 1) - ^(3) - ^(1) 



(3.53) 



(3.54) 



Note that the expression in the square bracket coincides with the leading twist anomalous 
dimension 77V (2.67) up to the color factor 2Cf- Proportionality || W jv| | 2 ~ 7^ is certainly 
not accidental. Exploring this interesting connection goes, however, beyond the tasks of 
this paper. 

3.4 Example: N = 1 

In the simplest case N = 1 one obtains 



(i) 



(2) 



9 \T/ V 1 — —A. 

Inserting these expressions in Eqs. (2.70) and (2.72) yields 

3ig 



(3) 
N=l 



(3.55) 



(dO) N=1 
Tn=i 



(nn) 
3ig_ 
(nn) - 



(Oa - Qi) + (Q2 - Qi) 
(Q 2 - Qi) - (Q2 - Qi] 



z k =0 
z k =0 



-6iq L gF +fi ^q L , 
-QqLgF+^^qL 



(3.56) 



where q L = Ul- 75)5. 
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The first equation in Eq. (3.56 is exactly the contribution of left-handed quarks to the 
identity (1.2) [21, 22] which is easy to derive by direct computation. The coefficient is 
correct because in our normalization 

In its turn, the operator Tjv=i appearing in the second line in Eq. (3.56) is familiar from 
studies of power corrections ~ l/Q 2 to the deep-inelastic scattering [39]. Its anomalous 
dimension equals jt ± = 4/3(iV c — 1/N C ) [39] and indeed coincides with the leading-twist 
anomalous dimension 7at=i (2.67), in agreement with our general result. 4 

4 Intermediate summary: kinematic projection operators 

For the reader's convenience we reiterate here the main result. 

As explained in the introduction, calculation of kinematic corrections ~ t / Q 2 , m 2 / Q 2 
to hard reactions in off-forward kinematics requires taking into account contributions of 
higher-twist operators that can be reduced to total derivatives of leading-twist opera- 
tors Oat. In particular picking up the contribution of the divergence of leading- twist oper- 
ators (dO) 7v is complicated because its matrix elements on free quarks vanish. Using QCD 
equations of motion this operator can be expressed in terms of quark-gluon operators as 
follows: 

2(nh)(dO) N = j^^[((4 N ,^))-((^ N ,Q))]+... , (4.1) 

where ~C§ = ^Q\,Q 2 -,Q^ and Q = {Qi,Q 2 ,Q3} are the light-ray operators with "good" 
conformal transformation properties, Eqs. (2.53) and (2.54), the scalar product ((. . .)) is 
defined in Eq. (2.69) and the coefficient pn m Eq. (2.40). Explicit expressions for the "wave 
functions" W tv are given in Eqs. (3.45), (3.46), (3.47). The ellipses stand for contributions 
of quasipartonic operators which are irrelevant in the present context. 

Using this representation and orthogonality of the "wave functions" of different oper- 
ators with respect to our scalar product we derive the following "kinematic projections" 
for the complete set of flavor-nonsinglet non-quasipartonic light-ray operators: 



i 9 Q{zi,z 2 ,z 3 ) = (nh) E \m \\2 ~*Nk(zt,z 2 , z 3 ) d\{dO) N + ... , (4.2) 



N=l k=0 

igQ(z 1 ,z 2 ,z 3 ) = (nn) V J2 PNk + I {-1) N+1 ~$ Nk (z 3 ,z 2 , Zl )d%{dO) N ^ 



where ||^at|| 2 is given by Eq. (3.54), the coefficients pNk are defined in (2.75), (2.48) and 
the functions Wjvfe(^i, z 2 , z 3 ) are given by 

~$ Nk ( Zl , z 2 , z 3 ) = (S+) k V ' N (zi, z 2 , z 3 ) , (4.3) 

where the differential operator S+ depends implicitly on the conformal spins jx,j 2 ,j 3 of the 
components of ~$ n, cf. Eq. (2.74). The ellipses in Eq. (4.2) stand for the contributions of 



4 It appears that the coincidence 7^ = "/n=i remained unnoticed over the 30 years. 
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other twist-four operators which are not reduced to total derivatives of leading twist ones, 
hence their matrix elements are not reduced to leading twist parton distributions. 

As a byproduct of this study, we have found a series of "genuine" twist-four flavor- 
nonsinglet quark-antiquark-gluon operators which have the same anomalous dimensions as 
the leading twist quark-antiquark operators: 



2{nn)T N 



(N + l) 



+ ... 



(4.4) 



These operators are not reduced to total derivatives, so that they contribute e.g. to the 
total cross section of deep inelastic scattering. 



5 T-product of electromagnetic currents (I): Twist expansion 

In the remaining Sections we apply the general framework developed above to the calcula- 
tion of the time-ordered product of two electromagnetic currents 

jfi(x) = q(x)j^Qq(x) , (5.1) 

where Q is a matrix of the quark electromagnetic charges Q = diag{e u , e^, . . .}. Our goal 
will be to obtain the complete expression for 

T ilv {z x ,z 2 ) = iT{j^( Zl x)j u (z 2 x)} , T afi ^ = a^a^T^ (5.2) 

to the twist-four accuracy, including the higher-twist operators related to total derivatives of 
the leading-twist ones. Here x^ is a four- vector and z\, z 2 are real numbers. Schematically 

z 2 ) = T*= 2 ( Zl ,z 2 ) + T*= 3 ( Zl ,z 2 ) + T l =\ Zl , z 2 ) + ... (5.3) 

Conservation of the electromagnetic current implies the Ward identities 

d»T^( Zl ,z 2 ) = Z2 \iP»,T„ v ( Zl , Z2 )] , 

d v T^( Zl ,z 2 ) = Zl \iP v ,T^( Zl , Z2 )] . (5.4) 

In addition, translation invariance along the line connecting the currents implies that 

e izP - x T^( Z1 , Z2 ) e~ izP - x = T^{ Z1 + Z , Z2 + Z ). (5.5) 

Both relations, Eqs. (5.4) and (5.5), are only valid in the sum of all twists but not for 
each twist separately. For the Ward identity, this property was noticed in Refs. [11-14]. 
Implementation of the electromagnetic gauge invariance beyond the leading twist accuracy 
has been at the center of many discussions. By contrast, violation of the translation 
invariance condition (5.5) for the contributions of a given twist has never been emphasized, 
to the best of our knowledge. As a manifestation of this problem, the structure of the twist- 
four contribution depends in a nontrivial way on the positions of the currents. For example 
we find that the kinematic twist-four part of the T-product with symmetric positions of the 
currents, iT{j^(x)j u (— x)} t=4 , is considerably more complicated compared to the case when 
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o 
o 
o 



a) 



b) 



Figure 1. Leading-order contributions to the time-ordered product of two electromagnetic currents 



one of the currents is taken at the origin, iT{j^(2x)j u (0)} t= . In the calculation of physical 
observables the difference must be compensated by different structure of the Nachtmann- 
type kinematic corrections that originate from the subtraction of traces in leading-twist 
operators. An example will be given in Sec. 8. The lesson is that the distinction between 
finite-i and target mass corrections oc t/Q 2 , oc m 2 /Q 2 that originate from the operators of 
different twist has no physical meaning. In particular the existing estimates of kinematic 
corrections to DVCS from the contributions of twist-two operators alone can be misleading. 

The leading-order contributions to the T-product are given by the two Feynman di- 
agrams shown in Fig. la,b (and crossing-symmetric diagrams which are not shown for 
brevity) : 

T$(z 1 ,z 2 ) + ... (5.6) 

The ellipses stands for quasipartonic operators of twist-four and operators of higher twist 
t > 4. 

A simple calculation yields (see e.g. Ref. [37]) 



2ir 2 



^^-q(zix)^^ v Q 2 q{z 2 x) + (/J, -H- v, Z\ •<-)■ z 2 \ 
c z 12 v ' 



Zl,Z 2 ) 



fl 



SlT 2 X 2 Zl 2 J 
+ ( fj, -H- V, Z\ O z 2 



dux a q{z l )Q 2 lll F a0 (z^) l5 + i{2u - l)F c 



i p ivq{z 2 ) 
(5.7) 



The Wilson lines between the quarks are always implied. Going over to spinor notation we 
obtain 



T ( tM,z 2 ) 

apap 



7T 2 X 4 zh 



{ X af3 ^a( z l x )^f}( z 2x) ~ Xp-(z 2 x)Xa(z\x) 

+ Xp-a\xa(zix)xp(z2x)-i)p(z2x)-ip a (zxx) | (5.8) 
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and 



1T 2 X 2 Z 12 J 



duj^Oi) u(S/)^ a (z2i) + «(a;/) a /j(^i) #( z 2) 

-Xa(^l) u(%f)ap(Z2l) +u(xf) / 3a(z% 1 ) Xp(z 2 ) 

+ X p(z 2 ) [u{xf)p a {z% 2 ) + u(xf) a$ (zf 2 )] x*{zi)} , (5.9) 

where for brevity we omitted the matrix of quark charges Q 2 . 

The operators in Eqs. (5.8), (5.9) do not have a definite twist yet. The twist separation 
is our first goal, where, in difference to the existing results (e.g. [37]) we have to take into 
account operators containing total derivatives. 

5.1 Handbag diagram, Fig. la 

For definiteness, let us consider the operator built of the left-handed quark and antiquark: 

O a a{ziX, Z 2 x) = ^ a (z\x) [z\X, Z 2 x]^ a (z 2 x) . (5.10) 

Our task in this section is to expand O aaL in contributions of different twist 

O a a(z 1 x,z 2 x) = 0^ 2 (zix,z 2 x) + O t =£(z 1 x,z 2 x) + O t ( ^(z 1 x,z 2 x) + ... 
To this end we start from the Taylor expansion 



(5.11) 



O aa (zix,z 2 x) 



CO ^ 

w, 



k=0 



2 k kV 



...x akak O 



aai...afc,Qcii...afc 



where 



O 



(k) 



Here and below 



Note that the operator ak Afc is symmetric under interchange of the pairs of 

indices (aj, &j) (ay, ay), i,j = l,...,k. 

In order to simplify the notations we will use multi-indices Ak^A^: 



(Z1,Z 2 ) = 3 ai&1 . . . d ak a k O a a{ziX, Z 2 x) 



x=0 



n 



d 

! dx^ 



aa : M 



(5.12) 



(5.13) 



(5.14) 



so that 



We also write 



A k = {a,ai . . .a k } 



A k = {a,a\... d k } , 



A k A k ~ ^aat...a k ,aut...a k \^i^) 



X a X ai ...A 



"A- 



9 



o 



8X A * d\ a 5A Q i ' ' ' dX a * 



(5.15) 
(5.16) 

(5.17) 



and similar in the "dotted" sector. 
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5.1.1 Leading-twist projection operator 

An important advantage of the spinor formalism for the twist separation is that sym- 
metrization and trace subtraction in Lorentz indices are replaced by one requirement: 
picking up the leading-twist contribution corresponds to symmetrization separately in dot- 
ted and undotted indices. The symmetric (alias leading-twist-two) part of (5.13) can be 
written as 



A,A=0 



where 



and 



O fc (A,A) = \ Ak X Ak 0f A ( Zl ,z 2 ) = (\d\) k O ++ (z l x,z 2 x) 



k^-k 



2=0 



++ (z\x, z 2 x) = X a X a O a a(zix,z 2 x) = ip + (zix)[zix, z 2 x]^ + (z 2 x) . 
Inserting (5.18) into (5.12) we obtain 



(5.18) 
(5.19) 
(5.20) 



Oifiz,, z 2 ) = d a d a X; [{k {dxdf [ [(XdX) k O ++ (z lX , z 2 x)\ x=q \ a _ =q . (5.21 



Here 



d d 
(dxd) = — / 



(XdX) = X a d aa X a 



dX« dX c 
Note that we are using shorthand notations 





dX a ' 



dX £ 



(5.22) 



(5.23) 



for the derivatives with respect to the auxiliary spinors, whereas d aa = ^-dn stands for 
the derivative over the four- vector x^. We hope that this similarity will not lead to a 
confusion. 

The sum in (5.21) can be rewritten as follows 



1 



Oi= 2 (z l5 z 2 ) = dJaJ2 [{k + lW (^ d ) k [ E ¥ ^^9XrO ++ (z 1 x,z 2 x) 



k=0 

oo 



m=0 



x=0 



d a da T7r ^—y I (dxd) k O ++ (z 1 n, z 2 n) 



A,A=0 



A,A=0 



(5.24) 



where n aa = X a X a and we took into account that only the terms with k = m survive. 
Finally, using that d aa (dxd) k+1 = 2(k + l)d a d a (dxd) k one obtains 



1 



O t aa 2 (ziX,Z 2 x) = -d a a / ^ E 



I oo 



(dxd) 



fe+i 



++ (zin, z 2 n) 



1 f 1 

-d aa Tl(x,X) / duO ++ (uzin,uz 2 n) 

2 Jo 



a,a=o 



(5.25) 
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where the projector A) is given by the following expression 

[Uf](x) = U(x, A)/(A, A) = f; ^£/(A, A) | (5.26) 

Note that the operator 0++(uzin,uz2n) on the r.h.s. of Eq. (5.25) "lives" on the light 
ray, n = A A, thus n 2 = 0. Hence the leading-twist part of the nonlocal quark-antiquark 
operator at x 2 =f=0 is expressed in terms of the light-ray operator. 
II(x, A) is indeed a projector since 

n 2 = n n(x, A)n(n = A ® A, rf) = n(x, ??) . (5.27) 

It satisfies a number of relations that are useful in practical calculations: 

n(x, A) = , 



n(x,A)A 7 3 Q = —-{xd) ia n(a3, A) , 



n(x,A)A 7 5o, = -(s9)^Q,n(a;, A) , 



1 - f 1 

n(x,A)A a A Q = n(x, A) x 2 B aa / duII(ux,A). (5.28) 

2 Jo 

These relations can easily be verified making use of Eq. (5.26). 

It is instructive to compare the expression for Tl(x, A) with the leading-twist projection 
operator in the conventional vector notation. It can be shown that (see Appendix C) 



U(x, A)/(n = A A) = £ P k (x, d y )f(y)\ y=Q , (5.29) 

fc=0 

where 



^) = i(-V) fc/2 C« f -#=1 (5.30) 



2 k k\" y ' k \^xW / 

and C^\x) is the Gegenbauer polynomial. This form is not unique. Equivalent represen- 
tations can be found in Refs. [15, 17, 37], see also the book [40]. The above expression for 
the projector II in the spinor formalism is much simpler and, therefore, easier to handle. 

5.1.2 Twist three 

The next task is to separate the twist-three contribution, 0^{z\x, z 2 x). To this end one 
has to contract one pair of indices in the operator (5.13) and symmetrize over all the others. 
Since 

OliOtjp) . q . _ . pp. 

c u a.i<Xi u OLjOtj ^0LiO.j U 

a contraction involving two derivatives gives rise to a twist-four operator. The only pos- 
sibility to get twist three is therefore to contract the quark field index with the one of a 
derivative. Hence one can construct two operators 

e afca [9 ai a! • • ■ d ak a k O a a{zix, z 2 x)]\ x=0 , t aak [d ai6ll ■ ■ ■ <9 Qfcfife O adl (zix, z 2 x)] \ x=0 , (5.31) 
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and symmetrize over all remaining open indices multiplying with the auxiliary spinors: 



O 



(A, A) = (XdX) k ~ 1 d a+ O k a + (z l x,z 2 x)\ x=0 



(M> ( 

o^" 3) (A,A) = (\d\) k ~ 1 d + aO k+ & (z 1 x,z 2 x)\ x=0 



(5.32) 



where for brevity we do not show the dependence of , on ^ ne coordinates z%. 

Since the operators in (5.32) are the only existing independent twist-three operators, 
the twist-three part of the local operator with open indices can be written as 



O 



(t=3) 
A k A k 



TT /*}(*-3) , -pr /r)(t-3) 

L1 A k A k u k ^ LL A k A k u k 



(5.33) 



Explicit expressions for the coefficients LT^ ^ , 11^ ^ can be found from the requirement 

(t=3) 

that O A £ is symmetric under the permutation of the pairs of indices (c^, on) (<x,-, 

z, j = 1, . . . , k and reduces to either 3 ^ or 3 "* by a contraction of one pair of indices 
and multiplication with auxiliary spinors for the remaining ones. One obtains 



JJ £)(*-3) 



p + l)!F 



n 



(t=3) 



P + l)!] 2 



ft, 



J'=l 



t=l, 



\ 



oi*" 3) (A,A) 



A=0 



J'=l 



t=l, 







(t-3) 



(A, A) 



A=0 



(5.34) 



It remains to rewrite this result in terms of nonlocal operators. Let 

0^= 3 \ Zl n,z 2 n) = [d a+ O a + )( Zl n,z 2 n), 
(i=3) (zm,z 2 n) = [d + aO + & ]( Zl n,z 2 n). 

Inserting (5.33) in (5.12) and using explicit expressions from (5.34) we end up with 



(5.35) 



O t c £{z 1 ,z 2 ) = -J2 



k(dxd) k 



k-l 



(dx) a daO it=3) ( Zl n,z 2 n) + {xd) 6l d a O { t =3) {z 1 n 1 z 2 n) 



(dxd) k 



E 



X a daO {t - 3) (z in , z 2 n) + XadaO^ (zw, z 2 n) 



A,A=0 



A,A=0 



(5.36) 



Finally, using 







Bx^ 



^(BxB) k {Bxf = 2(k + 2)(BxB) k B 6 







(dxd) k (xdji = 2{k + 2){dxd) k 8 a 



(5.37) 



this expression can be rewritten as 



O t a 3 (z 1 x,z 2 x) 



du uu 



un) 



da 7 U(x, A)A 7 A a O (t " 3) ( Zl un, z 2 
+ 9 a ^II(x, A)A^AQ,d) ( * =3) (2inn, z 2 un 



(5.38) 
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where IL(x, A) is the projection operator defined in Eq. (5.26). We remind that n aa = X a X a 
and application of the projection operator involves sending A, A — > in the final expression. 

As the next step, we have to unravel the twist-three light-ray operators on the r.h.s. 
of Eq. (5.38) to separate the contributions of interest. Omitting EOM terms that contain 
D a+ ip a or D + aip a one obtains 

O i - t=3) (zin,z 2 n) = zi[D a+ i) + \(zin)ip a (z 2 n) - igzi2 \ du z^ip+izin) f ++ (z 21 n) ip + (z 2 n) , 



o 



£> ( * 3) (zin,z 2 n) = z 2 ^ a (zin)[D + aip+](z 2 n) - igz u \ du z 21 ij) + (zin) f ++ (z 21 n) ip+(z 2 n) , 

(5.39) 



o 



where the terms in /++, /++ originate from differentiation of the Wilson line, cf. Eq. (3.12b). 
Such terms give rise to quasipartonic twist-three operators which are not relevant in the 
present context, so that we will discard them altogether. The remaining contributions (i.e. 
without the integral) can be rewritten in terms of the momentum operator (3.3), e.g. 

^X)0^= 3 \z in ,z 2 n) = z 1 (^A)[iP Q+ ,^+(z 1 n)V' a (^n)] 

= z 1 \iP ++ ,ip + (z 1 n)ip_(z 2 n)] - z 1 [iP-+, ip + (z 1 n)ip + (z 2 n)] .(5.40) 

The second term on the r.h.s. of (5.40) has the desired form: it is a total transverse 
derivative of the leading- twist operator. The first term can be related to descendants of 
twist-two operators using the following identities: 

d 1 

H a -Q^O ++ {zin,z 2 n) = + ^{zin,z 2 n) + -Zi[D^ + i) + \{zin)^ + {z 2 n) 

+z 2 d Z2 + -(zi,z 2 ) + ... , 
\iP- + ,0 ++ (z 1 n, z 2 n)] = [D- + il) + }{z 1 n)^ + (z 2 n) + 2d Z2 + -(z 1 n, z 2 n) + . . . , (5.41) 

where 

+ ~(zin, z 2 n) = $ + (zin)'ip-(z 2 n) . (5-42) 

and the ellipses stand for the contributions of quasipartonic operators. Excluding the 

term [D \-ip + ](zin)ip + (z 2 n) from Eqs. (5.41) one obtains an equation for the operator 

+ -(zin, z 2 n) 

d Z2 z 2 iO + -(zin,z 2 n) = (fid\)0 + (zin, z 2 n) - y [iP_+, + (zin, z 2 nj\ (5.43) 
which is easily solved to give 

(zin, z 2 n) =(fid\) / duO + {z\n, z^n) / du[iP |_, + (z\n, z^n)] . (5.44) 

Jo 2 J 

Using this expression in (5.40) we obtain 

0i f Zl r 

A Q C ( * =3) (zin,z 2 n) = -zi[iP a+ , ++ (z 1 , z 2 )] H / dw iP ++ ,d a O ++ (z 1 ,w) 

z 2 f zi r r n 
--^/ dw iP a+ ,0 ++ ( Zl ,w) (5.45) 

2^12 ho L L J J 
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and similarly 

A A 0( i=3 )(zm,z 2 n) 



~ 7 2 



-Z 2 \i~P+a,0 ++ (z 1 ,Z 2 )] + — / 

z 12 Jz2 



z 2 

2z 12 







/ dw 




> Z2 





iP ++ ,9aO ++ (u),Z2) 

?;P +tV ,0 ++ («;,Z2)ll • (5.46) 



Finally, inserting these expressions in (5.38) one obtains after some algebra 

_ If 1 f Zl dw r 

O t c ^ i (zix,Z2x) = -3 Q q,I1(x,A) / duu — iP ++ , ziO ++ (ziu,wu) + z 2 ++ (wu, z 2 u) 



Z2 



1 



+ -9„ 7 n(x,A) / duu 2 



Z\ 2 




r 


dw 


1 z 2 


Z\2 


r 


dw 




Z\2 



i(Pn) Q7 , ziO ++ (z\u, wu) 



(5.47) 



Note that the operators O++ on the r.h.s. of this expression are on the light cone, i.e. 
++ (wu, z 2 u) = ++ (wun, z 2 un). Making use of Eqs. (5.28) one can move all A-dependent 
factors to the left of the leading-twist projector Tl(x, A) and rewrite the answer in terms of 
the operator 



O l + +(zix,z 2 x) = Tl(x, \)0 ++ (zin, z 2 n) = [UO ++ ] (z±x, z 2 s 



(5.48) 



i.e. the leading-twist part of the nonlocal quark-antiquark operator ip + {zix)tjjj r {z 2 x) off 
the light cone. In this way one obtains 



O t a ^ 3 (z 1 x,z 2 x) = - I duu I — 



dw 



212 



iPfj_, (xa ll d) a c t zi0 1 + +(ziu,wu) + (xa^ d)o ia z 2 O t + +{wu, z 2 u 



t=2, 

++ 1 



+ In ud a aX 2 d u 



iP u , z\0\l{z\u, wu) + z 2 O t + +(wu, z 2 u) 



t=2. 



+ . . .(5.49) 



which is our final result. It is easy to check that 

x aA C£= 3 (zix,z 2 :E) = d ai *O i =£(z 1 x,Z2x) = 0. 

Note that the terms in the last line of Eq. (5.49) are by themselves twist-four. They 
must be included in order to achieve the proper separation of twist-three and twist-four 
contributions (i.e. their role is to subtract unwanted twist-four contributions that are 
present in the other terms), but can be omitted if the calculation is done to the twist- 
three accuracy only. The ellipses stand for the contributions of twist-three quasipartonic 
operators and EOM terms. 



5.1.3 Twist four 



Taking into account that Qfe & (5.13) is symmetric under permutations of the 

pairs of indices (aj, dj) (ay, ay), i,j = 1, . . . , k, one has three possibilities to project a 
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twist-four operator: 



pAkAk _ \ e aia2 e ai&2 (ja ya.z A afe )(A"A <i3 \ ak ) 
pA k A k _ e aia gdid (\ a 2 A Qfc )(A" 2 X° k ) 

pA k A k = eai a e a 2 a( X a 2 _ _ _ ^(A^A" 3 . . . A° fc ) , (5.50) 



so that 



\x=0 ' 



p 1 ■ OW = -(Xd\) k - 2 d 2 ++ (z 1 x,z 2 s 

p 2 .o( fc ) = -{xd~x) k - 1 d a aO a '"(z 1 x,z 2 x)}\ 

P 3 ■ = -(XdX) k ~ 2 d a+ d +A O aA (z 1 x,z 2 x)]\ x=0 . (5.51) 

However, because of the Fierz identity 

d 2 ++ = {XdX)d a& O a& - B +a d +6i O a& (5.52) 

only two projections are independent, e.g. P3 • can be eliminated in favor of the two 
others. Introducing the operators 

j<m i^{j,m} 
k 

j'=i 

we can write the twist-four part of the quark-antiquark operator with free indices as 

{t =f = Y n ( , a) , M ab P b • . (5.54) 

A k A k A k A k 
a,b=l 



The matrix of coefficients M ab can be found applying the projection operators p^ kAk ; 
P 2 Afcj4fe to th 
lation gives 



pA k A k ^ o e q Ua ^j on anc [ re q U i r i n g that the result reproduces Eq. (5.51). A short calcu- 



(5 - 55) 

Inserting this representation, Eq. (5.54), in Eq. (5.12) and using that 

j ^ 9 « 9n (i) = --jfe x 2 ^^{dxd) k ~ 1 , JJ s d « a «n( 2 ) = -it ^^(SxS)*- 1 (5.56) 

one obtains after a little algebra 

U2 



1 7 1 



A: 

oo 



A,A=0 



~2. 



A,A=0 

(5.57) 
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where 



0^ A \x;z 1 ,z 2 ) = [d fj$ OW](z 1 x,z 2 x), 
0^(x;z u z 2 ) = \ a \«[d 2 x O a< i}(z 1 x,z 2 x) 



(5.58) 



The result in Eq. (5.57) can be rewritten in an integral form which is more convenient in 
applications: 



O 



z 1 ,z 2 ) + 0^(z 1) Z2), 



(5.59) 



where 

{ a H(zi,z 2 ) 

( a m (zi,z 2 ) 



tin 



du 



(1 + lnu) x a a H — lnux 2 d, 



(x;uzi,uz 2 ) 



\nux a a + -(^—^ — l+\au 



{x; uz\, uz 2 ) , (5.60) 



Using Eq. (3.12a) the operators (j4) ,0 (B) can be expressed in terms of quark-gluon oper- 
ators. Neglecting the quasipartonic contributions one obtains 



OW{n- Zl ,z 2 ) 



~i 



dw w 



F + F 



(zi,w,z 2 ) , 



(s) (n; zt,z 2 ) = z 1 z 2 \iP l _ l ,[iP^,0 ++ (z 1 ,z 2 )]] -2z 12 z 2 F(zi, z 2 , z 2 ) + z x F(z u z u z 2 ) 



+ 



dwlyZi(w - z 2 ) V(zi,w,z 2 ) - 2d Zl [F + F](zi,w, z 2 ) 
+ z 2 {w - zx) V{zx,w,z 2 ) - 2d Z2 [F + F](zi,w,z 2 ) 



}, (5.61) 



where F, F, T>,T> are defined in Eq. (B.l). Using Eq. (B.13) the result can further be 
rewritten in terms of operators of the conformal basis (2.53), (2.54). 

5.2 Gluon emission diagram, Fig. lb 

The gluon emission diagram shown in Fig lb is written in terms of quark-antiquark-gluon 
operators, Eq. (5.9), and does not contain twist-two terms. The leading twist-three con- 
tributions to this diagram are due to quasipartonic operators of the type which 
are irrelevant for the present study, as we repeatedly stressed above. The goal is, therefore, 
to isolate the subleading twist-four contributions. The expression in Eq. (5.9) contains 
eight terms which are all similar so that it is enough to work out twist separation on one 
example. 

Expanding the nonlocal operator in the Taylor series 



oo ^ 

i)a{zix)f la {wx)ll)p(z 2 x) = ^ — 



n=0 



2 n n\ 



i=l 



Ipa (Ziy)fya {wy)lpp {z 2 y) 



y=0 



the task is reduced to picking up the twist-four part of the local operators 



- 0/3701... «„,ai...a„ 



n 9 <xi«i 



i=i 



^a{z\y)f ia {wy)iljp{z 2 y) 



(5.62) 



(5.63) 
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By construction, T Q( 3 7ai ... Q , li « 1 ...a n is symmetric inaf>7 and under the interchange of the 
pairs of indices (c^di) (ctjCtj). Hence there exist three apriory inequivalent possibilities 
to contract one pair of indices and symmetrize over the other ones: 5 



Pi ■ T — e l3a X y X a A Ql A a "T Q ,^ 7Q , 1 . i . Qni Q, 1 ... Qn , 
k=l 

n 

P 2 .T = e ^A^A a A d A dl J] X a ^T aPiai ... anM ... an , 

k=2 
n 

P 3 -T = e^A^A^A" 1 J] \ ai \^T aMai ... anAl .. Mn . (5.64) 

fc=2 

Obviously, the contraction of two "undotted" indices as in Pi, P2 and the contraction of 
two "dotted" indices, P3, give rise to tensors with different Lorentz spin. It is easy to check 
that P3 results in operators that do not contribute to (dO) and can be omitted. Hence 
only the first two projections, Pi and P2 are relevant. 

The rest of the calculation is very similar to the separation of higher-twist contributions 
to the quark- antiquark operator. Let 

n 

III = (e,s a <9 7 + da JJ d aj B a . , 

n n 

n 2 = ^epajdad-fdadaj Y[d ai B ai . (5.65) 

Note that the operators have the same symmetry properties under permutation of in- 
dices as T. Suppressing the open indices, the twist-four part of the tensor T Q!/ 37»i...a n ,ai...a n 
can be written as (cf. Eq. (5.54)) 

t<«> = ' K 1 M ■ (5.00) 



(n + 3)[(n + l)!] 2 V A ' V \ -2 3 J^-T^ 

It is easy to check that Pi ■ [T — T^ =4 -*] = i.e. T^* -4 ** is indeed the twist-four part of the 
tensor T. 

Inserting (5.66) in (5.62) and taking into account the identities 



xff[x**°*IL 2 = "^ry (x^(xd)d aa - \x 2 d^d aa ^ (dxd) n+l (5.67) 



3 The e aai contraction can be rewritten in terms of Pi, P 2 because e mi X = e alS X ai + e 0ai X° 
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one finds 



x^[ipa(z 1 x)f ia (wx)ilj l3 (z2x)] = - -x^d^I^x, A) J du F(uz 1 , uw , uz 2 ) 

+ £ (Xapdpa + Xppda&J A) J du (1 - U 2 ) 

x ^z\uV{ziu, wu, Z2u) + 2{wd w + 2)F{uz\,uw, WZ2) j . 

(5.68) 

Similarly, for the operator that involves the anti-selfdual gluon field one obtains 



1 



x a ' y [ilJc t (zix)f^^(wx)^i3(z2x)} = --x a pdpoJl(x, A) J duF{uz 1 ,uw,uz 2 ) 

+ \ ( X a/3 d P°< + Xaadp^ U(x, A) ^ du (1 - U 2 ) 

x iz 2 uD(z\u, wu, Z2u) + 2{wd w + 2)F{uzi,uw, j . 

(5.69) 

The operators T,J-,T),T> are defined in Eq. (B.l). 

6 T-product of electromagnetic currents (II): Kinematic twist-four 

Adding together the contributions of the two diagrams in Fig. 1 we can write the twist-four 
contribution to the T-product of two electromagnetic currents in the following form: 

2 



apap 



2™4^3 



12 



A(x; zi,z 2 ) - A(x; z 2 , z\) 



x a odpaE(x; zi,z 2 ) - XfSad a Mx; z 2 , z{) 



x 2 x /3 pd a aC(x;z 1 ,z 2 ) - x a ad jj pC{x;z2,z 1 ) } 



+ 



(6.1) 



where the ellipses stand for the contributions of the right-handed quarks and quasipartonic 
(four-particle) operators. Note that this expression is explicitly symmetric under the inter- 
change of the pairs of spinor indices (a, a) (/3$) and, simultaneously, the replacement 
z\ <H> z 2 . 

The invariant functions A, B and C are given by the leading-twist projection of nonlocal 
quark-antiquark-gluon operators and can be written as 



A(x; zi,z 2 ) = II(cc, X)A(n; Z\,z 2 ) 



(6.2) 



(and similar for the other two), where A{n; z\, Z2), B{n; z\, Z2), C(n;zi,z 2 ) are light-ray 
operators, i.e. they only involve the fields on the light-cone. The ^4-term originates from 
the handbag diagram in Fig. la only, the C-term is due to the contribution of the gluon 
emission from the hard propagator, Fig. lb, and the 5-term receives both contributions. 

In order to get explicit expressions, we have to collect the results of Sec. 5.1.3 and 
Sec. 5.2 and rewrite them using operators from the conformal basis, Eq. (B.13). This 
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calculation is described in Appendix D. At this stage it proves to be advantageous to use 
the relations in Sec. 2.5 to rewrite all contributions of Q\ and Q3 operators in terms of Q 2 . 
Remarkably, all results can be written in terms of the light-ray operators R, R that involve 
one and the same integral over the light-cone coordinate of the gluon field: 



R(zi,z 2 ) 
R(z 1 ,z 2 ) 



(nh) 
(nh) 



Z2 
Zl 



dw (w - z 2 )Q 2 (zi,w,z 2 ) 
dw{z\ -w) Q2(zi,w,z 2 ) 



(6.3) 



This property is only true for the sum of both contributions in Fig. 1 but not for each 
diagram separately. 
We obtain: 



A(n;zi,z 2 ) 



B(n;z\,z 2 ) 



C(n;zi,z 2 ) 



i J du ju 2 \q.uz\z 2 [iP 1 * [iPjj, 0++]] (z\u, z 2 u) 
+ f z 2 d Z2 — \nuz 2 d 2 z Z\ 2 ) R(uzi,uz 2 ) 



Z\2 
~2 



zid zi - - lnu zid Zl z 2 i j R(uzi,uz 2 ) > , 

^ / ^| lu 2 (l-u 2 +u 2 lnu) ziz 2 [iPV\iP fl ,0 ++ }}(z 1 u,Z2u) 
o Jo u I 



(1 - u ) z 2 d Z2 + (1— « +tt Inn) z 2 d l z z\ 2 



+ 



212 

(1 - u 2 ) ( zid Zl — I + (l-u 2 +u z lnu) Zid zi z 2 \ 

V z 2l 

du 



R{uz\,uz 2 ) 
R(uzi,uz 2 )\, 



[R(uzi,uz 2 ) + R(uz 2 ,uzi)] 

o In U 



(6.4) 



The last step is to obtain a "kinematic projection" for the operators R(z\,z 2 ), R(zi,z 2 ). 
Since the contributions of Q± and Q3 operators are no more present, we only need to use 
the expansion (4.2) for Q 2 which involves the functions 



^1(Z U Z2,Z 3 ) = (S^ U) ) k ^(z 1 ,Z2,Z 3 ) 



(6.5) 



(1,3,1) 

We remind that S + ' 2 ' 2 is the three-particle step-up operator for conformal spins j\ = 

(2) 

l,j2 = 3/2, j'3 = 1/2. Explicit expression for ^ v N J (zi, z 2 , £3) is given in Eq. (3.46). The 
possibility to represent the result in a simple form is based on the following identity (see 
Appendix E): 



dw{w-z 2 )(sf U) f*%\z l ,w,z 2 ) = {S^' 0) ) k [ Zl dw(w-z 2 )^ ) (z l ,w,Z2), (6.6) 



where S*.,. ' 2 ' 2 = z 2 d Zl +2zi+w 2 d w +3w+z 2 d Z2 +z 2 and S\_' ' = z 2 d Zl +2zi+z 2 d Z2 . In other 
words, a simplification arises because the quark-antiquark-gluon operator gets integrated 
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over the light-cone coordinate of the gluon field in the expression for the T-product of two 
electromagnetic currents with a particular weight, Eq. (6.3). 

Making use of the explicit expression for Eq. (3.46), we get 

dw(w-z 2 )^ ) (z 1 ,w,z 2 ) = r N (z 1 -z 2 ) N+l , (6.7) 

where 

r N = ~ \(N + 2) [iP(N + 3) + il,(N + 1) - ^(3) - ^(1)] . (6.8) 
Substituting this expression in Eq. (4.2) one obtains (cf. Eq. (2.49)) 

R{Zl,Z 2 )= 2^}_^ — .. ||2 [S\' ') z 12 + d + {do) N + ... 

N=lk=0 ^ 7V " N >> 

OO OO 

= - E E ^(^' 0) ) fc 4 +1 d%{dO) N + . . . , (6.9) 

N=l fc=0 

and similarly 

OO OO 

R( Zl ,z 2 ) = - £ E 9 k + (dO) N + ..., (6.10) 

N=l k=0 

where the ellipses stand for the "dynamical" contributions of "genuine" quark-antiquark- 
gluon operators that we will omit hereafter. The coefficients uJNk are defined in Eq. (2.50). 

It remains to substitute the expressions (6.9), (6.10) into Eqs. (6.4) and perform the 
integration over u, which becomes trivial. 

Let 

*%r ) (z 1 , Z2 ) = (S^) k z». (6.11) 

Using the identities 

b a,± Z 12 - z 12<\± > d z\± - b o,± d z i b - i2 ) 

the final results can be written in terms of two functions, Nk (z\, z 2 ) and ^^ k 2 (z\,z 2 ), 
where the first one is already familiar from the expansion of twist-two light-ray operators, 
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^\z u z 2 ) = y N - k 2 \z 1 ,z 2 ), cf. Eq. (2.49). We obtain 



B(n;z 1 ,z 2 ) 



AA(n;z 1 ,z 2 ) = A(n,z 1 ,z 2 ) - A(n,z 2 ,z 1 ) 
1 f 1 

= ~.z\z 2 / duu 2 lnit [iP M [iPfj., ++ (z±u, z 2 u) - ++ (z 2 u, z\u)]\ 
4 Jo 

7V,oddfc=0 V A ' y 

+ -^^2, 1 z 2 (a 21 + 5 22 )< 1) (z ll , 2 )} a* (ao)^ , 
1 r 1 

-z\z 2 I du(l -v 2 + it 2 lnit) [iP 11 [iP IJ ,,0 ++ ]] (21 it, 2 2 i0 
° Jo 

1 ^ r 3 1 31 

- 4 JCE (Jy+a) (^ (y+t+ a) {"'iff'^.*) - ^'(o.*.) 

-ff^[ 2 ' +Z2+ i( JV+fc+4 )^< 2 '*.- 22 ' 9 «)]*™ , ^' 22 )}^< 80 »«' 

.00 

cWi^) = t^i 2 £ E7FTWTiy^ )(zi ' Z2) ^ (90 ^ ' (6 - 13) 

TV.odd fc=o v """A 

which is the final result. The functions ^ Nk (z±, z 2 ) and ^ {z\, z 2 ) can be written more 
explicitly with the help of the integral representation in Eq. (2.52). In this way the results 
can be brought to the form similar to Eq. (2.26), cf. [31]. 

6.1 Light-ray operator representation for R, R 

Kinematic projection for the operators R(z\,z 2 ) (6.9) and R(z\,z 2 ) (6.10) and as a con- 
sequence the results in Eq. (6.13) are written as Wilson OPE in contributions of local 
operators of increasing spin and dimension. In applications to hard exclusive reactions a 
nonlocal representation in terms of light-ray operators is often preferable, as it allows one 
to calculate physical amplitudes directly in terms of parton distributions rather than their 
moments. For this purpose we need to find a light-ray representation for R, R operators 
which we can rewrite as 

00 00 j 
= -zi 2 Y,T,WT^ { m\^^)d k + {dO) N , 

N=l k=0 

oo oo 

R(zuz 2 ) = -z 12 Y^Y.WT^ { Nk\^^)d k + {dO) N . (6.14) 

N=l k=0 

As the first step, using \iP &a , <9 Q <9 d O n ] = (N + l) 2 (dO) N and applying [tP + , - ] k to both 
sides one obtains after a short calculation 

\iP &a ,d a B & 8 k + O n ] = {N + \) 2 d k + (dO) N + h(2N + k + 3)[iP dQ! , \iP a& , O^On}} . (6.15) 
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Multiplying both sides by coNk^^\zi, Z2) and summing over N, k this becomes 
^u Nk (N + l)H Nk ( Zl ,z 2 )d k + (dO) N = 

Nk 

= [tP &a ,d a d & ++ (z 1 ,z 2 )] - ^S^ 1 \iP dttX ,[tP a a,,0 ++ {z 1 ,z 2 )]], (6.16) 

where we took into account that k(2N + k + 3)wNk = ^Nk-i and used Eq. (2.49) to rewrite 
the sum over O+On in terms of the light-ray operator ++ {z\, z 2 ). 

The expression on the l.h.s. of Eq. (6.16) can be transformed into the sum for the 
R(z\,z 2 ) operator, Eq. (6.14), using the following identity: 

The simplest way to verify this relation is to observe that the integral operator 

[K<p](z u z2) = C da f dp t <p(z? 2 , 4) . (6.18) 
Jo Jo P 

intertwines the representations T jl=1 <g> T j2=1 and T jl=3 / 2 ® T j2=1 / 2 , that is 6 

KS^ = S { }^K. (6.19) 
As a consequence, since ^^)(z\,z 2 ) = (S^) k z^ 2 one derives 

(6.20) 

which is nothing but Eq. (6.17). 

Applying this operator to both sides of Eq. (6.16) we obtain 

R(z 1: z 2 ) = -z 12 J\a £ d(3^\iP^,d^O ++ (zf 2 ,z^)] 

- i4 fi) [^ 7 J^77.0++(% a 2.4i)]]} (6-21) 

and, similarly 

R(z u z 2 ) = -z 12 C dp f 'da"([iP^,a 7 ^O ++ (0f 2) 4))] 
Jo Jo a I 



l 5 (_M) ti p77 5 [ i p 7 .,0 ++ (^ 2 ,z 2 3 1 ))]]| . (6.22) 



The relations in Eqs. (6.21), (6.21) are between the light-ray operators, i.e. they involve 
fields "living" on the light-ray n 2 = 0. The OPE of the product of two electromagnetic 



6 A general approach for the construction of two-particle intertwining operators is described in Ap- 
pendix B in Ref. [30]. 
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currents (6.1) involves the leading-twist projector II acting on the light-ray operators (6.2). 
For the calculation of matrix elements it is convenient to rewrite 



U(x;zi,z 2 ) = U(x,X)R(zi,z 2 ) , TZ(x;zi,z 2 ) = n(x, X)R(zi, z 2 ) 



.23) 



directly in terms of 0++{z\x, z 2 x) = Tl(x, X)0 ++ (zin, z 2 n) (5.48) i.e. in terms of the 
leading-twist part of the nonlocal quark-antiquark operator i tp + (zix)ip+(z2x) off the light- 
cone. Using 



1 



Il(x, \)d,ydjO ++ (zi, z 2 ) = -d^{zid Zl + z 2 d Z2 + l)0% _l_ {z\x, z 2 



(6.24) 



we obtain 



K(x; Zl ,z 2 ) = z 12 I da f° dd\ S^^, [iP„, Of 2 (z? 2 x, z^x)]] 



Jo 



[3 



(3 



^2 2' _ 1 > 



P 



a I (——) 

"J cA2 2>\ 



K{x;z 1 ,z 2 ) = z l2 \ dp J da-^Sl^'liP^liP^O^iz^x^^x)}} 



Substituting these expressions in Eqs. (6.4), (6.2) gives the kinematic part of the twist- 
four contribution to the T-product of two electromagnetic currents (6.1) in terms of the 
leading-twist operator 0++(ziX, z 2 x). To save space, let us introduce the notations 



1 {z 1 ,z 2 ) = [iP^,[iP^,0^_ 2 (zix,z 2 x)]] 
2 {z u z 2 ) = \iP^,-^O t = 2 {z 1 x,z 2 x) 



(6.26) 



and 



A0 1 (z 1 ,z 2 ) = 1 (z 1 ,z 2 )-0 1 (z 2 ,z 1 ) 
A short calculation gives: 



A0 2 ( Zl ,z 2 ) = 2 (z u z 2 ) - 2 (z 2 , Zl ) . (6.27) 



AA = — J du < u 2 In u z\z 2 AO\{uz\, uz 2 ) — u 2 f da f dfi z 2 2 AO\(uzf 2 , uz 2l ) 
4 Jo I Jo Jo 



I daf 



o Jo 



zl^^AO^uz^uz^ + z^S^AOxiuz^uz^) 



i 

u z In u I da 



o 



z 2 S$* ) AOi {uzf 2 , uz 2 ) + z x S^ 2 } Ad {uzi , uz% x ) 



u[l - S(u)] f da r% \z^A0 2 {uzf 2 , uz^) + z{ 2 A0 2 {uz{ 2 ,u Z2 \\ 
Jo Jo P 
-1 

u(l + \nu) I da z 2 A0 2 (uzf 2 , uz 2 ) + z\A0 2 {uz\, uz 21 ] 
o L 



(6.28) 
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1 



™1 pot 

- I dui (1— u 2 +u 2 hiu)ziZ20i(uzx, uz 2 ) + (1 — u 2 ) I da dj3 zf 2 Oi(uzi 2 , ^21) 
I Jo Jo 

+ (l-u 2 ) I da r^lz^S^O^uz^uz^ + z^S^Oxiuz^uz^ 
Jo Jo P 

- (l-u 2 +u 2 lnu) [ da \z 2 S^ ) 1 (uz% 2 ,uz 2 ) + z 1 S ( }^ ) 1 {uz 1 ,uz% 1 
Jo L 



,/9 



-(1+li 2 ) / da I i- Z^O^UZ^U^) + -212^2(^12,^21/ 



'0 JO /3 

H — (1 — u 2 h\u) I da z 2 2 (uzf 2 , uz 2 ) + z\0 2 (uzi, ^21) 



(6.29) 



C = -^12 J duj^ da j^dplis^ 'AO^uzf^u^) - (<5(iZ) + ^)a0 2 ( 
where 

AA(x; zi, z 2 ) = A(x; zi, z 2 ) - A(z; £2,^1) • 
To arrive at these expressions we used the following identities: 

rl pot ^ pi pot 

d Z2 z i2 da d(3-= 0(zi 2 x, z 2l x) = ~ da d(3 0(zi 2 x, z 2l x) 
Jo Jo P Jo Jo 



a 13 \ 
UZ±2, UZ 2 i) 



(6.30) 
(6.31) 



d 2 Z2 z 2 2 [ daTddoizf^z^x) 
Jo Jo P 



daO(zf 2 x, z 2 x) 



(6.32) 



which can be verified by examining their properties under conformal transformations. 

Our final expression for the T-product of electromagnetic currents to the twist-four 
accuracy is given in the next Section. 



7 T-product of electromagnetic currents (III): Final expressions 

In most of the discussion so far we considered the contribution of left-handed quarks ip only. 
Adding the contribution of the right-handed quarks, x, is straightforward and reduces to 
a redefinition of the leading- twist operator O++. From now on it will be implied that 



with 



OX 2 { Zl ,z 2 ) = Tl(x,\)0 ++ (z 1 ,z 2 ) 



0++(z 1 ,z 2 ) = il) + {z 1 n)'ip + {z 2 n) - x+(z 2 n)x+(zin) 
1 



= 2 [l( z l n )i0- - l5)q(z 2 n) - q{z 2 n)iji(\ + 75)^1™)] 
= -[O v (z 1: z 2 ) - O v (z 2 ,Zi) - O a (zi,z 2 ) - A (z 2 ,zi)] . 
Oy and Oa are the usual vector and axial-vector light-ray operators 

O v (zi,z 2 ) = q{zin)i/hq(z 2 n) , O a (zi,z 2 ) = q(zin)ft'y 5 q(z 2 n) 



(7.1) 



(7.2) 



(7.3) 
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Note that Oy always enters antisymmetrized and Oa symmetrized over the quark positions. 

Now we are in a position to write down the final result. The T-product of two elec- 
tromagnetic currents including all contributions of leading-twist quark-antiquark operators 
and their total derivatives to twist-four accuracy can be cast in the form 



L Q/3d/3 



(zx,z 2 ) 



TT 2 X 4 Z 



12 I 



+ X 



X 8 gd a aC(z 1 ,Z 2 ) ~ X a ad g gC(z 2 ,Zi) 



+ 



(7.4) 



or, equivalently, in vector notation [31] 7 



(IV 



TT 2 X 4 zf 2 I 



ie 



+ x- 



(x^dv + Xvd^X + (xpdv - x u d^)Y j , 

(7.5) 



where 



^/i(zi,z 2 ) = ©^(zi,^) -®/i( z 2)^i) + x IJi AA(zi,z 2 ) 
W„( Z1 ,Z 2 ) = z 2 ) - <B^ 2 , z x ) , 



X(zi,z 2 ) 
Y(zi,z 2 ) 



(zi, z 2 ) - C(z 2 , zi) 
(zi, z 2 ) - C(z 2 , zi) 



(7.6) 



%$pa(zi,z 2 ) = a^^B^zi, z 2 ), S^ avf 3 = g m g v p + g va 9^ ~ 9piu9ap, and a totally antisym- 
metric tensor is defined such that €0123 = 1- 

The invariant functions AA(z±, z 2 ) and C(zi,z 2 ) are given by the expressions (6.28) 
and (6.30), respectively, with the redefined operator O++ as shown in Eq. (7.2). Both 
these contributions are twist-four. In addition, the function ^S a a{z\,z 2 ) contains all twists 
starting from the leading one: 



BodOzi, z 2 ) = K~*{zi,z 2 ) + W-£{z u Z 2 ) + %^{Z U Z2) + ■■■ 
Collecting all expressions, we obtain: 



(7.7) 



If 1 

®Ld {Z1,Z 2 ) = -d a a 



du 0\+ (uz\x, uz 2 x) , 



W^{zi,z 2 ) 



1 



— / udu I — 

Jz 2 Z\ 2 



dv 



iP M , (xa^d) a aZ 1 O t ^(ziu,vu) + (xa' 1 d)o l aZ 2 OW(vu,z 2 u 



+-lnud a aX 2 dpp 



ipW, zxO^Hzm, vu) + Z2OX+ (vu, z 2 u) 



®Ld (*1>*2) = X 2 d a aM( Zl ,Z 2 ) 



(7.8) 



with M(zi, z 2 ) given by Eq. (6.29). We have checked (using the representation in Eq. (6.13) 
for the twist-four functions AA, B, C) that this expression satisfies the Ward identities (5.4) 
up to twist-five terms. 

7 The function W 3 in (7.5) corresponds to in the notation of Ref. [31]. 
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As mentioned in Sec. 5, translation invariance relation (5.5) is only recovered in the 
sum of all twists as well. As a consequence, dependence of the results on the positions of 
the electromagnetic currents, Z\ and z 2 > is nontrivial. We have found that the expressions 
for higher-twist contributions can be simplified considerably for a special choice, when one 
of the electromagnetic currents is at the origin, e.g. z\ = 1 and z 2 = 0. In this case we 
obtain 



®jf*(l,0) 



and 8 



<B*= 3 (0,1) 



2^= 4 (0, 1) 



du I dv 







J ppva ~ ^ppvu 



x p (f + In u d^x d v 



-x <9 M / du I dv 



3u 2 , / uv 
1 + — + -ln (1--- 

u v V uv 



Oi(«,u) 



+ 



n 



?1 



1+- 3+~+- + 



U V 



1 + v 2 



lull-™ 

U V 



2 (u,v)\ (7.9) 



^= 2 (0,1) = J\uO^(0,u) 



du / dv 







x p (f + In n <9 M a; d v 



iP",0^(u,u) 



-x I du I dv 



3u 2 , / uv 
1 + — + -ln (1--- 



Oi(u,u) 



+ 



u 



u 



U V 



1+- 3+~+- + 



U V 



1 + v 2 



lull-™ 

U V 



2 (v,u)}. (7.10) 



Also 



AA(1,0) = -/ du dv{ --AOi(«,u) + - + In 1 - -- ))A0 2 (u,v) } (7.11) 



U V 



u V 



and finally 

C(1,0) 



fin / dv 







+ 



i> 2u 2 / it i> 
_ In 1 

V u V \ uv 



AOi(u,v) 



V u 



V u 



V 1 u 



___ 2 + - + -- --In I--3 



v 2u 



1 



u v 



U V 



C(0,1) 



du / cfo< - AOi(v,u) 



u 



u 
u 



1 + 



In 
2n 



A0 2 («,u)l, 



(7.12) 



where AOi t2 (u,v) = Oi >2 (u,v) - Oi >2 (v,u). 

8 The results for 25 M (1,0) and 23 M (0, 1) are related by hermitian conjugation, ©^(21,22) = (03^(22, zi))^, 
taking into account that 0^(u,v) — (0++(v, u))* and 01,2(11, v) = (£'1,2(1', u)) t . 
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Operator product expansion in terms of nonlocal light-ray operators Oi^{v,u) has an 
advantage that the matrix elements are given directly by parton distributions rather than 
their moments. The above expressions still appear to be rather complicated, however. 
A possible drawback of this representation is also that the operator O2 = iP 11 , d^O^r 2 
contains both terms of the type ~ d C^...^ and ~ d >l O flfll ,,,p n , cf. (6.15), (6.16), and one 
may like to have them separated. 

This can be achieved by using another representation which is usually referred to as 
conformal OPE (see e.g. the review [27]), in terms of integrals over light-ray positions of 
local conformal operators, similar to Eq. (2.26). To this end we define the (axial)vector 
conformal operators as 



0T\y)=(d Zl+ d Z2 ) N cf f'|^ 



V ( A )(zix + y,z 2 x + y) 



Zi=0 



(7.13) 



and 



(dO)T\y) = J^^^^iy)] = \i^,0^ m (y)]x^... X ^. (7.14) 



The twist-two contributions can be written as 



t=2 



asf 2 (i,o)-!8f 2 (o,i) = <)„yj 



7V,odd 



(N + 2)\ J 



t=2 



-^= 2 (l,0)-^= 2 (0,l)=dp J2 



2(2iV + 3) f 1 
(N + 2)\ 



duu N u N+2 



O v N {ux) 
Ol(ux) , (7.15) 



N,even 

where for arbitrary function the leading twist projection [. . .]u is defined as 

[f(x)] lt = [Uf](x) = U(x,X)f(n). 
The conformal expansion for the twist-three contributions reads 



(7.16) 



t=3 



2iV + 3 



7V,odd 



(JV + 2)!(JV + 1) 



2N + 3 



N,even 



(N + 2)\(N + 1) J 



+ dpx 2 d u 



N,odd 



2N + 3 f 1 

W+WJo 



duu N+2 u N I dvv iy+L lnv 



du{uu) N+1 iP u ,Ol{ux) 
du(uu) N+1 \iP u ,0$(ux) 

iP u ,0%(uvx) 



N+l 







2N + 3 



TV, even 



(N + 2)\(N + 1) J 



du{uu) N+1 iP u ,0^{ux) 



27V + 3 



ATodd 



(N + 2)\(N + 1) 



du{uu) N+l iP v ,Ol{ux) 



+ dpx 2 d u ^2 



2N + 3 



N,evcn 



(N + 2) 



duu N+2 u N dvv N+1 lnv 



iP",Ofr(uvx) , (7.17) 
lit 
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and, finally, for the twist-four terms one obtains in this representation a remarkably simple 
result [31] 

v m =4 = E {N + 2W + 2) Sl du (uU)N+1 { Xfl K^Mk 

+ ijV(JV + 3) J\v v N ~ l x 2 d„ [{dO) v N {uvx)] u , | , 



x " 4 = ^ (T+ + 2)i(jy + 2) jf* 1 (tifi) " (u ~ fi) I » 

Yt=4 = -^ ^ f/+ + 2)!(^+ 2 1 ) J/u^^+^J^dvv^ [{dO) v N {uvx)\ lt .(7.18) 



This representation reveals that the operator d 2 0^ 1 ,,,^ n which corresponds to [iP^^P^, Oat] 
in our notation, does not contribute to the answer for our special choice of the correlation 
function T{j p (x)j u (0)}. Thus, for this choice, kinematic power corrections to twist-four 
accuracy are related to the contribution of the divergence of the leading twist operators 
alone (apart from Nachtmann-type corrections). The T-product with symmetric positions 
of the currents, T{j p (x)j u (— x)}, includes both operators. The corresponding expression 
turns out to be much more cumbersome. 

8 Sample applications 

As we already mentioned, contributions of different twist in the OPE of the product of 
currents in off-forward kinematics are intertwined by electromagnetic gauge and Lorentz 
invariance. Implementation of the electromagnetic gauge invariance in off-forward processes 
beyond the leading twist accuracy has been at the center of many discussions, see e.g. [11]. 
By contrast, importance of the translation invariance condition has never been emphasized, 
to the best of our knowledge: the distinction between the kinematic corrections that orig- 
inate from contributions of the leading- and higher-twist operators is not invariant under 
translations along the line connecting the currents and does not have physical meaning. 
We want to illustrate this statement and also demonstrate restoration of the translation 
invariance in the sum of all twists on the simplest example: distribution amplitude (DA) 
of a longitudinally polarized p-meson. Our discussion closely follows Ref. [26]. 

To this end we consider a nonlocal quark-antiquark operator sandwiched between vac- 
uum and the /5-meson state [26] 

(0\u(z lX )^d(z 2 x)\p-( P )) = f p m p (e ■ x) / due'**™™ <f>«(u) + —z 2 l2 x 2 ^{u)+0{x A ) 

Jo 1 J-o 

(8.1) 

where m p and f p are the /9-meson mass and decay constant, respectively, p^ is the meson 
momentum, p 2 = m 2 and e p is the polarization vector, e • p = 0. A Wilson line is implied 
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between the quarks. The dependence on quark coordinates z\ and Z2 is fixed uniquely by 
the translation invariance condition 

{0\u((z 1 + 5)x)^d((z 2 + 5)x)\p~(p)) = e- l6( - px \0\u( Zl x)^d(z2x)\p- (p)) . (8.2) 

It is easy to verify that this relation holds for the above parametrization of the matrix 
element. 

The function <^>| (u) has a physical meaning of the momentum fraction distribution 
of the quark and antiquark in the longitudinally polarized p-meson and is dubbed the 
leading- twist DA. It is usually written as an expansion over Gegenbauer polynomials 

0||(ti l /x) = 6u(l-«) «Jv(/i)Cj /2 (2n- 1) (8.3) 

Af,even 

and is normalized by the condition du<p»(u) = 1 , so that ao = 1. The coefficients ajy in 
the Gegenbauer expansion correspond to the reduced matrix elements of local conformal 
operators On defined in Eq. (7.13). 

In turn, the function $>(u) is called the twist-four two-particle DA. Our task is to find 
all "kinematic" contributions to $>(u) which are proportional to m 2 [26]. For simplicity of 
the argument we inserted "^" between the quarks, in which case twist-three contributions 
drop out 

u{z 1 x)^d(z 2 x) = [u(z 1 x)^d(z 2 x)] t=2 + [u(z 1 x)^d(z 2 x)] t= ' 1 + 0{x 4 ) , (8.4) 

and, more importantly, will assume that (f>\\(u) has its asymptotic form at large scales, 
4>\\{u) = Qu{l — u), i.e. put all coefficients cin due to the conformal operators On with 
N > 2 to zero. In this case the matrix elements of twist-four operators (0O)n obviously 
vanish as well and we are left with a simple expression 

1 f 1 

[u(zix)gtd(z 2 x)] t ~ 4 = -x 2 zi z 2 \ dvv 2 [iP M , [iP^, [u(vz\x)$d(yz-2x)\i\ + ... (8.5) 
or, for the matrix element 

(0| [u(zix)ftd(z2x)] t=4: \p~ (p)) = — m 2 x 2 / dvv 2 (0\u(vzix)ftd(vz2x)\p~ (p)} + . . . 

4 Jo 

= -^/X(e • x)x 2 [ dvv 2 ! due- ivz ^ x U\\(u) + . . . (8.6) 
4 Jo Jo 

where the ellipses stand for dynamic twist-four contributions which are of no interest to us 
at present. 

One sees that the twist-four contribution alone does not have the form required by 
translational invariance. Moreover, it vanishes if Z\ = or Z2 = 0, i.e. if the quark of the 
antiquark field is at the origin. The invariance must be recovered in the sum with the term 
~ m 2 x 2 from the leading-twist contribution 

[u{zix)$d{z2x)] t=2 = [u(zix)ftd(z2x)]it = n(x, X)[u(zin)^d(z2n)] , (8-7) 
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where II is the leading-twist projector, e.g. in the vector representation (5.29) 

k=0 

oo 

= £ M M) 



n=0 



^ ^_2 (2; .a n ,) fc -V^ + o(x 4 )l/(n) 



fc=0 



71=0 



Here d n = d/dn p . One obtains after a simple algebra (cf. [26, 41]) 

[n(e • x)e- ipn ](x) = [(e • x) e - ipx ] Zi = (e • s) [e - * 1 * + ^pV f dvv 2 e - ivpx + 0( 



5.9) 



and therefore 



(0|[n(zix)^d(z 2 2;)] |p ) = f p m p (e-x) 



du 



-izf 2 (px) 



"°'<l>\\( u ) 



+ -mix 1 [ dvv 2 e- ivz ^ px \z't 2 ) 2 <l)\\{u) + 0(x 4 ) 



.(8.10) 



The second term (0(x 2 )) in this equation is the Nachtmann-type mass correction that 
originates from subtraction of traces in the definition of the leading-twist operator. Adding 
(8.6) and (8.10) one gets 



{Q\u{z\x)$d(z2x)\p ) = f p m p (e-x) / du 



i r 1 

+ -m 2 x 2 / dvv 2 e- i( - px)vz ^[(z^ 2 ) 2 - Zl z 2 ](p\\(u) + 0(x 4 ) 
4 ./n 



.11) 



Finally, using explicit expression <fr\i(u) = 6uu and rewriting [(^i 2 ) 2 — ^1^2] = z\2\z\2uu + 
(1 — 2u)z" 2 ], the second line can be simplified using 



du I dvv 2 e-^ px > z ^[(z^ 2 ) 2 - Zl z 2 ]uu = ^ / (uu) 2 e'^K 



.12) 



Thus, translation invariance is indeed restored and we obtain, comparing (8.11) with the 
ansatz in Eq. (8.1) 



$(«) = 12m 2 uV, 



5.13) 



which agrees with [26]. 

We emphasize that the the translational invariance is only rescued in our calculation 
for the asymptotic distribution amplitude <p\\{ u ) = 6uu, which is consistent with neglecting 
twist-four contribitions due to divergence of the leading twist (conformal) operators. For 
a generic <p\\{u) translation invariance will be recovered by adding the terms in (dO)]\f. 
The lesson is that kinematic corrections that originate from the contributions of operators 
of different twist in the OPE are all intertwined by Lorentz (and eventually also gauge) 
invariance. Estimates based on the contributions of leading twist alone can be misleading. 
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In the rest of this section we consider matrix elements of the operators in question 
sandwiched between proton states with different momenta, which are relevant for the de- 
scription of hard exclusive scattering reactions in Bjorken limit, such that the DVCS. To 
leading-twist accuracy, the relevant nonperturbative input is encoded in the generalized 
parton distributions (GPDs) which are defined by the matrix elements of nonlocal quark- 
antiquark light-ray operators. For example, for the vector operator (see, e.g. Ref. [8]) 



(p'\O v {z l n,z 2 n)\p) = / due- iP + {{z2 ^ +u)+Zl ^- u)) 



i 



x u{p') 

where 



ia +u A 

7+i?(u,£,t) + — z— ^(u,£,t) 
2m 



u(p), (8.14) 



P- T Af, A-^-p, t-A', (8.15) 

2 P+ + p' + 

m is the nucleon mass and u(p) is the nucleon spinor. As always, P+ = P^n^ etc. The 
GPDs H(u, £, t) and E(u, £, t) depend on Bjorken variable which we call here u in or- 
der to distinguish from the quark coordinate, skewedness parameter £ and the invariant 
momentum transfer t. 

Gegenbauer moments of the GPDs can be written in terms of reduced matrix elements 
of local conformal operators 

v 



j ducf H( U ,a,t) = f2 (-i) k m k - N [H Nk (t) + E Nk (t)} , 

fc,even 

I duC 3 J 2 (-) E(u,£,t) = (-l) k+1 m k - N E Nk (t) (8.16) 

fc,even 



where 



and 



{-i) N (p'\O v N ( n )\p) = u(p')Mp)HN(n;p,p') + < p > {p) E N (n;p,p') (8.17) 



N 



H N (n;p,p') = H Nk (t)A k + P?~ k , 

A;, even 
N+l 

E N (n;p,p') = ]T E Nk (t)A k + P? +1 - k ■ (8.18) 



fc.even 



Note that we define the generalized form factors Ej^^^) as coefficients of the scalar struc- 
ture u(p')u(p) instead of a more usual u{p')ia +v A v u(p) which leads to a redefinition of 
the H^k{t) form factors; hence a "tilde" in the notation. The description in terms of the 
generalized form factors of conformal operators has an advantage compared to the stan- 
dard parametrization in terms of simple moments, A n f,{t), B n i~(t), C n {t) [8], in that E^ k {t) 
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and H^ k (t) have autonomous scale dependence to one-loop accuracy. This property is 
convenient in applications [42, 43]. 

The distributions (8.14) are defined as matrix elements of the operators directly on 
the light cone, and they enter the OPE through the leading twist projection of the quark- 
antiquark operators at generic (non-light-like) separations: 

(-i) N {p'\[O v N {x)] lt \p) = u( P ')Mp) Hn(x; P , P ') + U{ - P ' )U ^ E N (x;p,p f ) ■ (8.19) 

m 

Using the explicit expression for the projection operator H(x,d n ) we obtain, to the 0{x 2 ) 
accuracy 



N 



N-2 



H N (x;p,p') = H Nk (t)(Ax) k (Pa 

fc,even 



4(iV + l) 



- y: (Ax) k (p X ) N - k - 



fc,even 



(k + l)(k + 2)tH Nk+2 (t) + {N- k)(N - k - 1) ( m 2 - - ) H Nk {t) 



N+l 



N-l 



E N (x;p,p') = Y E Nk (t)(Ax) k (Pa 



A:, even 



4(iV + 



- Y (Ax) fc (P, 



sAT-fc-l 



fc,cvcn 



+ l)(k + 2)tE Nk+2 (t) + (N + 1 - k)(N - k) ( m 2 - - ) E Nk {t) 



+ 2(N - k)m 2 H Nk (t) 



5.20) 



The terms 0(x 2 ) in these equations give rise to a finite-t and target mass corrections to 
hard exclusive scattering reactions, which is analogous to the Nachtmann correction in 
deep-inelastic scattering. For the particular case of DVCS such corrections were studied 
in [13-19]. The new contribution of this work is the calculation of corrections due to the 
"kinematic" twist-four operator (dO)^ which must be taken into account alongside the 
Nachtmann correction to maintain Lorentz and electromagnetic gauge invariance. 

The relevant matrix elements can be defined by an expression similar to Eq.(8.19): 



{-if-^p'\{dO) v N (x)\p) = u(p')Mp) H^(x) + ^u(p')u(p) E%\x) . 



1 



,(4), 



Taking into account that 



(p'\(dO)l(x)\p) 



1 

iV + l' 



d 
dxt 1 



{p'\O v N (x)\p) 



5.21) 



(8.22) 
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we obtain from Eq. (8.20) 



1 



N-l 




2(N + l) 2 



£(Ax)*(Px) 



N-k-l 



(k + l){2N + 2-k)tH Nk+1 (t) 



A;, odd 




> 



HNk-l(t) 




- ^2(Ax) k (Px) N - k (k + l)(2N + 2-k)tE Nk+1 (t) 



fc,odd 




These expressions have similar level of complexity compared to Nachtmann-type corrections 
in Eq. (8.20), but a somewhat different structure, which means that they may result in a 
different dependence on the quark momentum fraction. 

Applications to concrete processes require dedicated studies and this task goes far 
beyond the scope of our paper. Main question that has to be addressed is whether QCD 
factorization is applicable for a given process for kinematic twist-four contributions. 

9 Summary and Conclusions 

Operator product expansion belongs to the most important tools in the analysis of hard 
reactions in QCD. Whereas OPE for forward scattering (i.e. neglecting operators that 
include total derivatives) was thoroughly studied in connection with deep-inelastic lepton- 
nucleon scattering, the general case of off-forward kinematics has never been analyzed 
systematically beyond the leading twist. A new phenomenon which is specific for off- 
forward kinematics is appearance of higher-twist operators that are total derivatives of the 
operators of the leading twist. In this work we show that the structure of such contributions 
to the leading order in strong coupling is governed by conformal invariance. We call these 
terms kinematic, as they do not involve new nonperturbative input compared to the leading 
twist. Our study has been fuelled by potential applications to hard exclusive scattering 
processes that involve generalized parton distributions, in particular DVCS, transition form 
factors of the type 7* — > Mj, and to calculations of weak decays of i?-mesons in the 
framework of QCD factorization or light-cone sum rules. 

As a principal result, we provide a set of projection operators that allow one to extract 
the "kinematic" part of an arbitrary navor-nonsinglet chiral-even twist-four operator in 
QCD. We also calculate, to the same accuracy, the kinematic contribution to the time- 
ordered product of two electromagnetic currents. This result makes possible a complete 
calculation of the finite-t and target mass corrections in two-photon reactions that are 
amenable to factorization to twist-four accuracy. We stress that the proof of validity of 
QCD factorization in a given process, at least for the kinematic contributions considered 
here, is a separate question which is beyond the scope of this paper. We hope that fac- 
torization of kinematic twist-four contributions can be shown at least for some selected 
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reactions, because their structure is intertwined with the leading twist and strongly con- 
strained by electromagnetic gauge and Lorentz invariance. Clarification of this issue is 
certainly the most important task for future studies. 

Apart from this topic, our study can be extended in several directions. In this work 
we have only considered chiral-even operators, i.e. ones built of a pair of left-handed or 
right-handed quarks. For some applications, e.g. U-meson decays, it may be necessary to 
include contributions of chiral-odd operators (i.e. made of one left-handed and one right- 
handed quark) as well. This extension should be straightforward. A more difficult question 
concerns kinematic projection of flavor-singlet and in particular pure gluon operators. Such 
contributions are needed e.g. for the calculation of finite-i and target mass corrections that 
involve gluon GPD. From a theory perspective, it would also be interesting to reproduce 
some of our results using a different technique. Our final formula for the kinematic contri- 
bution to the T-product of two electromagnetic currents appears to be considerably simpler 
as compared to the intermediate expressions. This suggests that a more direct derivation of 
this result might be possible by considering three-point correlation functions of conformal 
operators with the electromagnetic currents to two-loop accuracy in the critical dimension. 
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Appendices 

A Evolution Hamiltonian for non-quasipartonic operators 

The matrix of anomalous dimensions for a set of local operators Oj is defined as 



where [0]i, Of are renormalized and bare operators, respectively. The renormalized oper- 
ators satisfy the matrix RG equation 



which, equivalently, can be cast in the form of an integro-differential equation for the 
generating functions (light-ray operators): 




[o]i = z£o£, 



(A.l) 




(A.2) 




(A.3) 



Here 7 is an integral operator, which we write as 




(A.4) 
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where a s = g 2 /Air. 

The operator Mqq is determined by one-loop counterterms to the nonlocal operator 
and can be written as 



M qq = N c M { q ) q - W 



1 



ir(-l) 



3C F - 2N c a q , 



(A.5) 



where a g = bo/4N c , &o = ^-N c — The "Hamiltonians" H^ 1 ) are 3x3 matrices which 
are written in terms of SL(2) invariant integral operators 



t(i) 



£ Q =(H 12 + H 23 )I + 



^12 ^23 ^ Z 12 



C^12 + ^12) 



e,(l) 
32 



-2ft 



12 



^ z^Uo -2(2nt 2 + nf 2 )J 



{2P 12 H'} 9 -H 



H13I + 



12ft 12 

5T7^r 2 no 



13 



2H7 2 + U d 1 , + P 23 uf 



1.03 



-2zi 2 ft 12 

^13^13 



^13 



(A.6) 



Explicit expressions for the two-particle kernels can be found in Ref. [30]. 

We want to find a scalar product such that Hqq is hermitian. Write, schematically 



: * 



#12 #13 

i^21 H 22 H 23 
yi^31 i?32 i?33 y 



(1 1 ^ 
(- - l) 



(A.7) 



where the superscripts on the three components of the "wave function" indicate the cor- 
responding SL{2) representations. Since the diagonal Hkk kernels are SL{2) invariant 
operators in their subsectors, it is natural to look for the scalar product of the form 



k=l 



where {&k, ^k) is the standard SU(1, 1) scalar product for the appropriate representation. 
We want to achieve that W = H, i.e. 



for arbitrary <£, \&. In particular we can take 









f 


\ 


$ = 




* = 




1) 




I / 




I 


/ 



(A.9) 



(A.10) 



and require that 



ai(0, #12^) (1,1,1) = a 2 {H 21 (f), i) 



(A.ll) 
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Since the evolution kernels are given by the sum of two-particle operators acting on a 
chosen parton pair, e.g. 

H 12 = Hg 2) + H{f ] + H[f } (A. 12) 

it is enough to consider two-particle kernels separately, and due to SL(2) invariance it is 
enough to consider their action on the corresponding eigenfunctions. For example, take 
H[f ] and (p = if> = (z 2 - z 3 ) n . Then 

H[? ] z2 3 = h X2 (n)z n 23 , flg 3 ^ = M*)*23 (A-13) 

and we have to require that 

01 h 12 (n) 1 1^23 1 1 (1,1) = a 2 h 2 i(n) H^ll^i) ( A - 14 ) 

which gives the relation between a\ and a 2 . Note that this expression has to hold for 
arbitrary n. Explicit expressions for hik can be read of Ref. [30] (see the results in invariant 
representation) after taking the color traces, or using the expressions in Eq. (A. 6). 

Completing this calculation for all cases one arrives at the result given in Eq. (2.69). 



B Derivation of Eq. (3.20) 

To save space, in this Appendix we will use the following notations: 

J 7 (z 1 ,z 2 ,z 3 ) = ig$ + (z 1 )f +a (z 2 )'ilj a (z 3 ,) , 

T(z 1 ,z 2 ,z 3 ) = igi) a (z 1 )fa + (z 2 )t/j + (z 3 ) , 

^(^1,^2,^3) = ig[Da+$+]{zi)f ++ (z 2 )if) a (z 3 ) , 

V(z u z 2 ,z 3 ) = ig^(z x )U + {z 2 )[D + ^ + ){z 3 ) . (B.l) 

We have to collect contributions to Eq. (3.11) from the operators containing two deriva- 
tives acting on the quark or the antiquark, and the contributions from quark-antiquark- 
gluon operators as indicated schematically by the second, third and fourth terms in Eq. (3.14). 
The contribution of the operator [D 2 ijj](zin)n'ijj(z 2 n) can be written as 

(zid Zl + z 2 d Z2 + 2)zi - -(d Zl + d Z2 )z\ [D 2 $\(zin)ni)(z 2 n) 

St 2 ~ [D 2 i)}{ Zl n)n^(z 2 n) . (B.2) 

As explained in the text, all terms in Sf 2 can be omitted since they vanish upon integra- 
tion (3.7). Adding the similar contribution with the derivatives acting on the quark field, 
one obtains for the sum of the second and the third terms in (3.14) 

h = -^d Z2 z 2 2 [D 2 if;}(zin)ni)(z 2 n) - \^d z ^ 2 ^{z\n)n\D 2 ^\{z 2 n) . (B.3) 
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Using EOM these operators can be rewritten in the form 



ip(zin)n[D 2 ip](z 2 n) = -2igip + (zin)f +a (z 2 n)?p a (z 2 n) = -2F(z 1 , z 2 , z 2 ) , 
[D 2 ij;](zin)ntj)(z 2 n) = 2igtp a (z 1 n)fa + (z 1 n)tf; + (z 2 n) = 2F{z x , z x , z 2 ) (B.4) 

so that one gets for (B.3) 

h = -d Z2 z 2 2 F(zi,zi,z 2 ) + d zl z 2 2 F(z 1 ,z 2 ,z 2 ) . (B.5) 

Quark-antiquark-gluon operators indicated schematically as ~ igipFip (the fourth term 
in Eq. (3.14)) actually come in two variations: with and without an extra covariant deriva- 
tive acting on the quarks. The contribution of operators without extra derivatives, call it 
I 2 , reads 



du 



(d ZlZl + d Z2 z 2 )z l2 - (d Zl +d Z2 )z 12 z 2 L 1 ^{zi)a»F fll/ {z% 1 )n u ^{z 2 ) 

du [ud Zl - ud Za ]zt a ${z 1 )ai>F ia/ (z$ 1 )n v il>(z2) , (B.6) 
where we can substitute 

${z 1 )v> l F lu ,(% 1 )n v 1>{z2) = -F( Zl ,z^,z 2 ) - F( Zl ,z^,z 2 ) . (B.7) 
The contribution of the operators with extra derivatives is (omitting the term in S± 2 ) 

h = - f du[ud zl z\ 2 ^{ Zl )nF^ + {z^)[D^^ 

(B.8) 

Rewriting the gluon strength tensor in terms of /, / one obtains 

$(z 1 )nF M .(z$ 1 )[D'*il>](z2) = -Mzi)(f +a (^)[D a + ^ + ](z 2 ) + f + ^ 1 )[D + ^ + ](z 2 )) 

= \f>{z u z^,z 2 ) - F^ ( Zl ,z l 2 \,z 2 ) - F^ (z 1} z^,z 2 ) (B.9) 



and similar 



1 



[I* i $](z 1 )nF M .(% 1 )il>(z2) = ^D{ Zl ,z^z 2 )-F^{ Zl ,z^z 2 ) 



F^( Zl ,z^,z 2 ), (B.10) 



where we used a notation F^ for the derivative over the corresponding light-cone coordi- 
nate: 



F^(z u z^z 2 )=d z F](z 1 ,z^ 1 ,z) 



Z=Z2 



F( 1 \z 1 ,z% 1 ,z 2 ) = d z F](z,z% 1 ,z 2 ) 



Z = Z \ 



and similar for F. 
Next, we write 



F {3) (zi,z^,z 2 ) = (d Z2 -—u—)F(zi,4 1 ,z 2 ) 
V zu du) 

FW( Zl ,z^,z 2 ) = (d Z2 - 2-u^\f{z 1 ,z^z 2 ) 



(B.ll) 
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and similar for T . Inserting (B.ll), (B.9), (B.10) in Eq. (B.8) and integrating by parts one 
arrives at the following expression for the sum I\ + I 2 + I 3 : 



(dO) ++ (zi,z 2 ) = d Zl zf 2 J 7 (zi,Z2,z 2 ) - d Z2 zf 2 J r (z 1 ,z 1 ,z 2 ) 



+ 



d Zl d Z2 z\ 2 



JF{z\,z\\,z 2 )+T{zx,z 2X ,z 2 ) 



d Z2 z\ 2 V{ Zll z^z 2 ) + d Zl zf 2 V( Zl , z^ 2 ,z 2 ) } + . .(B.12 



The ellipses stand for EOM, contributions of quasipartonic operators and terms propor- 
tional to Sf 2 which do not contribute to the projection (3.7) for the conformal operator. 

The last step, we have to rewrite the answer in terms of the operators from the con- 
formal basis (2.53), (2.54): 



(/_tA)(A/i)J r = ig[Q 2 (zi,z 2 ,z 3 ) - Qi(z 1 ,z 2 ,z 3 )] , 
(fi\)(Xfi)T> = 2ig[Q i 2 \z 1 ,z 2 ,z 3 ) - Q 3 (z 1 ,z 2 ,z 3 )] 
(fj,X)(Xp,)f = ig[Q 2 (z!,z 2 ,z 3 ) - Qi{z 1 ,z 2 ,z 3 )] , 

,z 2 ,z 3 ) - Q 3 (z 1 ,z 2 ,z 3 )] 



(B.13) 



where, as above, the superscript in , Q^p stands for the derivative in the corresponding 
argument. Using these expressions, rewriting the derivatives in Q 2 > Qp hi the form (B.ll) 
and integrating by parts we obtain the final answer given in the text, Eq. (3.20). 



C Leading twist projector 



Here we derive Eq. (5.29) which establishes a relation between the projectors in spinor and 
vector representations. By definition, the action of the projector (5.26) on the function 
f(x^) is given by the expression 



p/i(^)=E 



fc=0 



(dxd) k . I 



f[n^ = -Xa^X 



A,A=0 



(C.l) 



By a direct calculation one easily finds 

(Bxd)f(n) J-(3xo»X)-^-J{n) = (xd n )F(n) + ±(Xa" xa^d^fin) = x»K^f(n) 



where 

= {nd)d^ + <9 M - i^i^d 2 , 
and the derivatives are taken with respect to n. Next 
1 



{dxdff -Xa^X 



A=0 



[xKYf{n] 



n = = fOn)[ r K 



n=0 



(C.2) 



(C.3) 



- 60 - 



where 



Kfj, = 2n^(nd) - n 2 d^ + 2n h 



(C.4) 



is nothing but the generator of special conformal transformations for a scalar field ip with 
scaling dimension dim[</3] = 1 (<p(x) — > x~ 2 ip{x / x 2 )) . Thus 



[ l -xK 



f{d n ) [^-) exp 
n=o V oa ' 



f(9„) 



-ax^Ku 
2 * 

1 



1 



n=0,a=0 



\da J 1 — (xn)a + a 2 x 2 n 2 /4 



n=0,a=0 



(C.5) 



Taking into account that the last factor in the above expression is the generating function 



for the Gegenbauer polynomials C k 

1 



(i) 



k=0 



(!) f (xn) 



one obtains 



1 - (xn)a + a 2 x 2 n 2 /4: ^ k Vv 7 ^ 2 /? 



2 2 2 

-a x n 



k/2 



(C.6) 



f{dn) ( -xK 



^f{d n ) {x 2 n 2 fl 2 C^ ( JgL 



n=0 



n=0 2 k • 

Finally, inserting (C.7) in Eq. (C.l) one ends up with the expression in Eq. (5.29): 

p/iw = EiiW /2c i 1)ni5) 



(C.7) 



k=0 



Vx 2 d 2 



f(n) 



n=0 



(C.8) 



D Derivation of Eq. (6.4) 

Collecting the expressions derived in Sec. 5.1.3 and Sec. 5.2 we obtain 

1 f 1 f Zl 
A(n;z\,z 2 ) = — / duu 2 {\ + Inn) / dw w [J 7 + J 7 ] (z\u, wu, z 2 u) 

+ duu 2 lnu|ziz 2 [«P M [iP/j,, 0++]] {z\u, z 2 u) 
- 2z\ 2 ( Z 2 F(z\u, z 2 u, z 2 u) + Z\ T(Z\U, ZiU, z 2 u) 



+ / dw^z\{w — z 2 ) uD{z\u,wu, z 2 u) — 2d Zl \F + T](ziu,wu, z 2 u) 



+ z 2 {w — z\) uD{z\u, wu, z 2 u) — 2d Z2 [J 7 + J 7 ] [z\ 



u, wu, z 2 u 



}}, (D.1 
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B(n;zi,z 2 ) 



du(l- u 2 + u 2 lnu)^ziz 2 [iP M [?P M ,0 ++ ]] (z x u,z 2 u) 

z 2 F(z\u, z 2 u, z 2 u) + z\ T{z\U, Z\u, z 2 u) 



2z 12 



- f du / dwl — 2u 2 Inuw [J 7 + J 7 ] (z\u, wu, z 2 u) 
8 Jo Jz 2 



+(1 - u 2 + u 2 hiu) 



z\(w — z 2 ) [uV — 2d Zl J r — 2d Zl J 7 ] (z\u, wu, z 2 u) 
+ z 2 (w — zi)[uV — 2d Z2 F — 2d Z2 T\ (z±u, wu, z 2 u 



+2u z 



(w - z 2 )F + (w — zi)F 



(uz\,uw, uz 2 ) 



(1 — u 2 ) (w — z 2 )wd w F + (w — z\)wd w F (uz±,uw, uz 2 ) 



--(l-u 2 )u 



z{D + z 2 V 



(uzi,uw,uz 2 ) > , 



(D.2) 



C(n; z\, z 2 ) = — [ du(l — u 2 ) [ dw (z 2 — w)l —z±u V(uz\,uw,uz 2 ) +V{uz 2 ,uw,uz\) 
8 Jo Jz 2 1 2 



+ (wd w + 2) [lF(uzi, uw, uz 2 ) + T{uz 2 ,uw, uz\)\ > . 



(D.3) 



The next step is to rewrite the results in terms of the conformal operator basis, Eq. (B.13), 
e.g. 

-l 



C(n;zi,z 2 ) 



1 9 



~i 



16(nn)7 du{l ~ u2) j Z2 dw ^ Z2 ~ W ) 



z\d Zl Q 2 - uziQ 3 + (wd w + 2)[Q 2 - Qi] (uzi,uw,uz 2 ) 



nd zl Q 2 - uz\Qi + (wd w + 2)[Q 2 - Q{\ (uz 2 ,uw,uz{) > . (D.4) 



+ 



This expression can be simplified further by observing that in the present context we are 
not interested in contributions that are related to descendants of geometric twist-three 
operators, cf. Sec. 2.5. Consider the following operator identities that are easy to derive 
by explicit calculation: 

^[iP^+,^+(zi)f++(z 2 )^+{z 3 )} = Q 3 {zi,z 2 ,z 3 ) + d Z3 Q 2 (zi,z 2 ,z 3 ) + d Z2 Qi(zi, z 2 , z 3 ) 

-[iP + -,ip + (zi)f ++ (z 2 )ip + (z3)} = Q?,{zi,z 2 ,z 3 ) + d Zl Q 2 (zi,z 2 ,z 3 ) + d Z2 Q 1 (z 1 ,z 2 ,z 3 ) 

(D.5) 

The expressions on the l.h.s. are total transverse derivatives of twist-three quasipartonic 
operators. Hence they have autonomous scale dependence and cannot get mixed with 
derivatives of twist-two operators that we are interested in. In other words, the matrix 
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elements of these operators can self-consistently be put to zero so that Eqs. (D.5) allows 
one to eliminate all occurrences of Q3 (Q3) in Eq. (D.4) in favor of <3i,2 (Qi,2)- In the 
present case an additional simplification occurs: the operators Q\,Q\ drop out from the 
resulting expression and the derivatives acting on Q2,Q2 get collected into a combination 
z\d zi + wd w + Z2d Z2 + 3 which can be replaced by ud u + 3. Integrating by parts in u one 
arrives at: 



C(n;zi,z 2 ) 



1 a 



[nn) 



du \ dw{w — Z2) Q2{uZl,UW,UZ2) + Q2(uZ2,UW,UZl) . (D.6) 



Using the notations in Eq. (6.3) this expression becomes the one in the last line in Eq. (6.4). 
In the similar way we obtain after some algebra 



4 A(n; z±, z%) = z\Z2 / duu 2 lnu iP M , [iP^, ++ (uz\, 11Z2)] 



+ 



+ 



{nn) 



{nh) 



du u z 



-1 



dw 



~2 



Z2Q1 - wQ 2 



(UZI, UW, UZ2, 



-1 



lnu(z2d Z2 )z2i / dw Q2 — Q1 (uzi,uw,uz2) 



du u z 



dw 



ziQi-wQ 2 (uzi,uw,uz 2 ) 



lnu(zid Zl )zi2 / dw Q2 — Q1 (uzi,uw,uz 2 ) > (D.7) 



and 



8 B(n; z\, Z2) = z\Z2 \ du (1—u 2 +u 2 In u) {P^, [iP^, 0++(uzi, UZ2) 



+ 



+ 



(nn) 



jg 

(nn) 



du\(l-u +u \nu)(z2d Z2 )z2i I dw [Qi - Q 2 ](uzi,uw,uz 2 ) 



(l-u 2 ) / dw [z 2 Qi - wQ 2 ](uzi, uw, uz 2 ) 

J z 2 



du l(l-u 2 +u 2 \uu)(zid Zl )zi2 I dw[Qi - Q 2 ](uzi,uw,uz 2 ) 



dw \z\Q\ - wQ 2 ](uzi,uw,uz2) }■ 



~2 



fDi 



In this case the operator Q\ does contribute, but can be dispensed off using the relation 
(which is equivalent to Eq. (2.78)) 



Z\ — n 

Qi(zi,w,z 2 ) = d Z2 zi2 I dn- ^ Q 2 (zi, n, z 2 ) 



(zi — w) z 

With this substitution, after some algebra one arrives at the expressions given in Eq. (6.4). 
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E Identities for three-particle wave functions ^Nk{ z i: z 2 , £3) 



The "wave functions" v£j\nfc(zi, 22, ^3) (4.3) of three-particle nonquasipartonic twist-four op- 
erators have a rather complicated algebraic structure. In practice one usually has to deals 
with integrals of these functions over the light-cone position of one of the fields (usually 
gluon) with simple polynomial weights which follow from the structure of Feynman dia- 
grams. Such integrals can be much simpler compared to the wave functions themselves, see 
Eq. (6.6) for an example. Although this equation can easily be verified by a direct calcu- 
lation, it makes sense to show that it is due to conformal invariance. Apart of being more 
elegant, this derivation allows one to understand why a simplification has been possible for 
this particular case, and why it does not happen for some other integrals. 
Consider 



>N 



(z 1 ,z 2 )=f duu^P(zi,Z2i,Z2) = \- I dw (w - z 2 )^ < ${zi,w,z 2 ) , (E.l) 

J z 12 Jzi 



(2) 

where ^ N \z\, w, z 2 ) transforms according to the tensor product of the representations 
TJ' 1=1 ®TJ' 2=3 / 2 (g>T^ =1 / 2 under SL{2) transformations (2.16). If we prove that ^ N (z 1 ,z 2 ) 
transforms according to T J1=2 ®T J ' 2 , i.e. the integral operator in (E.l) is an intertwining 
operator between these representations, then 



(S^ 1) ) k $ N (z 1 ,z 2 ) = I duu (S^' 3 ' 3 ')**^ (zi,z^,z 2 ) (E.2) 
Jo L J 

which is equivalent to Eq. (6.6) since zf 2 (S+'^ ) k = zf 2 {S^ ) k . 

Behavior under conformal transformations can conveniently be analyzed using the 
reproducing kernel /C. Using the definition in Eq.(2.35) we can write 

*n(z 1 ,z 2 )= f duu III nPefc] (1,| '^^i(^i.fi)^|(^i»6)^(^^3) *!?(6,6,6)- 

fc=i 

(E.3) 

Now the integral over u can be taken, producing a product of two reproducing kernels with 
different spins, cf. Eq. (3.51) 

Jo 2 z 2 

which is nothing but a Feynman representation for the product of reproducing kernels, 
fully analogous to the product of propagators. 
Thus 

$ N (z 1 ,z 2 ) = \ [[[ n[^] (1,f,i) ^(^'6)^1(^1,6)^1 (zi,e 2 )/Ci (*2,6) 
fc=i 

(E.5) 

Following the analogy between reproducing kernels and propagators [44, 45], this equation 
can be represented as a "Feynman diagram" shown in Fig. 2, where the numbers in the cir- 
cles specify the conformal spins. Using SL(2) transformation properties of the reproducing 
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Figure 2. Diagrammatic representation of Eq. (E.5) 



kernels (3.29) it is easy to show that transformation properties in z-vertices are determined 
by the sum of the conformal spins of the propagators (reproducing kernels) attached to the 
vertex, in our case T J=2 for z\ and T J=1 for z 2 , which is the stated result. 

Thus, it is seen that the possibility to get a simple representation is tightened to the 
possibility to make use of the Feynman parametrization of the type in Eq. (3.51). Repeating 
the above arguments, one can easily derive 



du u 

du u 
l 

du 



o 



(S^ U) ) k ^ ) ](z 1 ,z^,z 2 ) 
{ S^) k ^}( Zl ,z^z 2 ) 



(2,lhk 



(S 



(2) 

duu^ K N '(z 1 ,Z2i,z 2 ) 



( 2 ' 2 )\fe 



duu^f 



(2) 

N (zi,z 21 ,z 2 ) , 



(S$' l] ) k r du^\ Zl ,z^,z 2 ), (E.6) 



however, e.g. 



duu 



(1,1,1) ^, Tf (i) 

N 



(zi,z 21 ,z 2 ) 



cannot be simplified in the similar manner. In our calculation of the OPE for the product of 
two electromagnetic currents this last integral appears in the contributions of the handbag 
and gluon emission diagrams, but cancels in their sum. This cancellation proves to be 
crucial in order to obtain a relatively simple final result, e.g. in the form (7.18). 
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